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The self-induced evolution of a vortex sheet was simulated by modelling 
the sheet using an integration of discrete elements of vorticity. 
Replacing small sections of a vortex sheet by flat panels of constant 
vorticity is found to reproduce more accurately the initial conditions for 
the Lagrangian simulation technique than replacement by point vortices. 
TLc flat panel method for the vortex sheet was then extended to model 
axisymme vortex sheets. The local and far field velocities induced 
by the asymmetric panels were obtained using matched asymptotic 
analysis, and some of the uncertainties involved in other models of the 
axisymmetric vortex sheet have been eliminated. One important result 
of this analysis is the determination of the proper choice of core size 
for a circular vortex filament which may replace a section of an 
axisymmetric vortex sheet Roll-up of both two-dimensional and 
axisymmetric vortex sheets was computed using the panel methods 
developed In the report. 
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Nomenclature 
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Roman Letters 

a Half-Width of a Panel (a « Vaw) 

b Semi-Span of the Wing or Radius of the Rotor 

(a 0 ,b 0 ) Position of Two-Dimensional Roll_,up Core 
E^ Elliptic Integral of the Second Kind 

F l Elliptic Integral of the First Kind 

Q n Panel-Position Polynomials; sec liquations (2.14) - (2.18) 

R Radius of the Axisymmetrie Vortex Element 

(R e ,c 0 ) Position of Axisymmetrie Roll-Up Core 
r Distance Between Points or to the Symmetry Axis 

r 0 Radius of Circular Vorticity Distribution 

r 1 Shortest Distance to the Vortex 

r 2 Longest Distance to the Vortex 

(r,z) Coordinates for an Axisymmetrie Cylindrical System 

s Line Integral Path Vector 

s Distance Along the Vortex Sheet or Panel 

1' Elapsed Time 

(u,v) Velocity Components for a Cartesian System 

v Velocity Vector 

w Width of the Panel 

X Q Distance from the Edge of the Sheet 

x Position Vector 

(x,y) Coordinates for a Cartesian System on the Vortex 


Greek Letters 

P Circulation 

T 0 Reference Circulation; Total Circulation 

e Curvature; Small Quantity Parameter 

6 Inclination of a Panel 

k "Gradient of Circulation"; Local Strength of the Vortex Sheet 

X Elliptic Integral Parameter; see Equation (2.24) 

to Vorticity Vector 

p Stream Function 


(Xl) Stream function Solution 
(Xe) Stream 1 'unction Solution 
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Chapter 1 
Introduction 

1.1. Foundation of Vortex Sheet Studies 

A common feature of flows at high Reynolds number is the formation of thin 

regions of sharp changes in velocity known as shear layers. In the limit of infinite 
Reynolds Number, the diffusive effect of viscosity is eliminated, and the regions of 
velocity gradients are reduced to surfaces of step changes in velocity. Such surfaces 
arc known as vortex sheets. Although the formation of shear layer is a viscous 
process, shear layers arc commonly approximated by vortex sheets in the study of 

invicid, irrotutional '’potential” flow. 

A vortex sheet can be represented as an integration over infinitesimal vortex lines. 

The dynamics of vortex litres, and thus the dynamics of the sheet, are given by 

Helmholtz’s vortex theorem. The theorem shows that vortex lines are material lines of 
the fluid and are convected by the local flow. In flows without sources or sinks, the 
distribution of vorticity determines the flow. If the flow is irrotational, the vorticity 
is concentrated in small regions of the fluid. Thus, an unsteady potential flow can be 
determined by tracking the vortex sheet as it is moved and deformed by the flow 
field. 

A rigorous approach to tracking a vortex sheet will involve solving singular integro- 
differential equations. The time integration can be simplified by using a numerical 
scheme involving finite increments in time. However even with this simplification, a 
rigorous representation of the vortex sheet will make the contour integral Intractable. 
Thus, a further numerical approximation must be used to represent the continuous 
vortex sheet using finite, discrete elements of circulation. 


The linear instability of a vortex sheet can be shown analytically, but its subsquent 
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development In difficult to descril*. Mxuutory experiments with thin (then* layers, 
which arc exacted to t« similar to vortex sheets, show spontaneous generation and 
growth of spatially periodic rollup, This instability of the shear layer, known ss the 
Kelvin Helmholtz Instability, was the subject of the first numerical study of vortex 
sheets by Roscnhcad [ij. In this study, a continuous two dimensional vortex sheet was 
replaced by a collection of discrete points of finite circulation, or point vortices. The 
initial flat geometry of the vortex ihect was represented by a Atrsight row of point 
vortices, Roscnhcad then adopted a Ugrnngian approach of following each point vortex 
over time. The motions of the vortices arc the result of the flow induced by the 
particular arrangement of vortices for each, step in time. Roscnhcad's calculations 
showed that the vortices form loose spirals within a small number of time steps if 
the initial arrangement of vortices were given a small periodic displacement. Because 
of the tedious manual computations involved, Roscnhcad proceeded no further than the 
early stages of rollup. As a result, the problems associated with this Ugrangian 
method were not known until many years later. 

The Lagrangian approach to the initial value problem involving vortex sheets is not a 
universally accepted concept. A pioneering effort by Birkhoff and Fisher [2] questioned 
the foundation of vortex sheet studies for invicid fluids. They were skeptical of the 
well posedness of applying invicid analysis to the dynamics of a viscously generated 
phenomenon. A recent review of the work in this area can be found in the article 
by Saffman and Baker [3]. 

Undaunted by questions regarding the wdl-posedness of the method, for the last 
fifteen years there has been an increasing interest hi refining Rosenheud’s technique. 
There are two reasons for the revival of interest in this area of study. First, the 

increase in availablity of powerful computers encouraged the development of 

computational fluid dynamics. Second, the desire to extend the invicid analysis for 

high speed aerodynamics to free vortex wakes has encouraged the development of the 

numerical approaches. 
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1.2. Instabilities Iri Nummieal Vortex Sheet Models 

The Jjagrpnglan approach to simulating * vortex sheet In often plagued by hiatanlltie* 
caused by the numerical approach, the manner in which vortex sheet* are modelled, 
and the inherent instability of an actual vortex sheet. The first la usually due to the 
inaccuracies inherent in numerical time integration techniques and In approximating the 
varying geometry of the vortex sheet. The second cant*; of Instability is due to loss 
of accuracy dictated by the amount of computation available for the problem. In 
order to satisfy the practical limitations in computing, the number and tlie complexity 
of vortex elements must be reduced. Since tlie I^grang^o technique uses the result of 
each time step as the input for the next time step, inaccuracies in the technique 
compound over time. However, indiscriminate elimination of instability In undesirable 
since the physical instability of a vortex sheet is indistinguishable from the artificial 
instability. The Inherent instability of tlie vortex sheet is the most important source 
of difficulty in the calculation of vortex sheet behavior. 

Vortex sheets can be modelled in many ways but whether these models can 
realistically represent an infinitesmally thin shear layer is often not clear. Moore [4] 
has demonstrated that point-vortex methods will always lie unstable. He found that 
the numerical Instability mimics the Kelvin Helmholtz instability in which the smallest 
wavelength resolvable, using points equal to twice the spacing, was found to be the 
most unstable. However, this instability causes the vortex sheet being represented by 
the point- vortices to cross itself, which is inconsistent with the allowed behavior of 
material surfaces. In the most successful point-vortex methods, the vortices are 

repeatedly redistributed along the sheet which unintentionally acts ns a low pass filter, 
damping out the instability to delay the onset of chaos. 

Baker [5] end Mnrman [6] have tackled the vortex flheet problems using a mixed 
Lagrangian-Eulcrian approach, commonly refered to as the "cloud-in-cdl” approach. In 
this scheme, the Euler Equation for the flow is solved for each cell formed by a 
spatial grid network. The use of the Euler Equation allows the use of a East Poisson 
Solver to significantly increases the efficiency of the numerical method. Murman 
confines the circulation to vortices which are convcctcd independently. Otherwise 
vorticity diffuses numerically, and structures with grid related length scales appear. 
Whether these structures represents physically realistic features which are too small to 
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Ikj revived by the grid In nut clew, On the other hand, although the email scale 
structures found l».y Moore ere suppressed, Murmau's ntethod tan not reproduce the 
smooth spiralling roll up expected In a vortex sheet. 

Spectral methods have been applied to a small numljcr of cases with dimple 
geometries to study vortex sheet stability. The problem addressed by Htwenhead has 
been studied using the Njjectrut method by Moore [7] and by Melron, llaker, and Uvs/ag 

lb], The former addressed infinite, smcl perturbations, while the latter addressed finite 
amplitude disturlmnecs. They both arrived at consistent result. Singularities In the 
form of infinite curvature of the vortex sheet were detected within a finite amount 
of evolution time. Although the method results appear to be extremely reliable, they 
arc too limited to be useful for practical problems such as the one discussed below. 
They do however illustrate the futitlity of rigorously following the evolution of a 
vortex sheet at infinltcsmal spatial scales. 

1.3. Application to wake of a lifting surface 

Any lifting surface of a finite span trails a shear layer In v.\Ac, commonly 
referred to as the wuke vortex sheet. This shear layer can be pnduccd by separation 
at the leading edge of a highly swept delta wing at high angles of attack or by the 
loss of vorticity from the vortex "bound” by the aerofoil. In the early 1920’s, Vraudtl 
hypothesized that the et fc es of wake vortex sheets roll-up into exponential spirals. 
Kndcn [9l using dimensional arguments, derived the behavior near the edges of the 
vortex sheet trailed by a wing with an elliptic lift distribution. This distribution is 
important in efficient airplane design since it induces a constant dowuwash which 
minimizes the induced drag for a given aspect ratio. Using dimensional arguments, 
Kaden found that the rollup produced a spiral of 2/3 power and from this the center 
of the spiral can be approximated. Betz [10] attempted to calculated the distribution 
vorticity in the spiral by conserving circulation and impulse. However, Betz’s constant 
which determines the exact shape and the distibution of vorticity with the spiral was 
later shown to be incorrect by Pull in [llj, 

Roscnhead’s method was first applied to "Kaden’s Problem" by Westwater [12]. In 
order to study the three dimensional geometry of the trailing vortex sheet, Westwater 
replaced the streamwisc coordinate by a time like coordinate to form an unsteady, two 


8 


dimensional process. This is equivalent to having the sheet crow-section being swept 
sway from the wing which generated it. This procedure, known as the Trefftz plane 
method, is now universally used to study wake vorlex sheets where the deformation 
in the streamwise direction is assumed to he amall. Westwater discretized the trailing 
vortex sheet by replacing 20 segments of the sheet with point vortices of the same 
circulation. The vortices are convectcd for each times tep by velocities they induce on 
each other. Although the vortices rolled up more or less smoothly, the spiral did not 
quite match Kaden’s result. The inaccuracy of Westwatcr's result was initially 
attributed to the insufficient number of vortices dictated by lack of computational 
devices. 

Attempts to improve Westwatcr’s results by using more point vortices have not been 
successful. Increasing the number of point- vortices was found to hasten the tendency 
toward chaotic displacement of vortices. Such failures in achieving higher resolution 
by increasing the number of point vortices arc attributed to the singular nature of 
point vortices. The velocity singularity at the .location of the point vortex can be 
removed by replacing an individual or a group of vortices with a nonsingular 
distribution of vorticity. However this method has been found to only delay the onset 
of numerical instability while introducing an ad hoc length parameter for the size of 
the distribution. The failure of the discrete vortex approach may be due to the 
impossibility of reproducing the sheet-like nature by a finite number of positions. On 
a sheet, each position posses its own unique orientation, which is the result of the 
sheet’s continuous nature. 

The method of equal-spaced rediscretization introduced by Fink and Soh [13, 14] is 
notable for its simplicity and for its success in delaying or eliminating instability. 
Fink and Soh found that replacing a vortex sheet segment with a point vortex at the 
middle of the segment will significantly reduce the discretization error. This in fact 
truncates the formula for the velocity associated with a vortex panel, a small segment 
of the vortex sheet. Ily placing point vortices at equal distances along the sheet, it is 
possible to replace equal length segments of vortex sheet by point vortices at the 
centers of the segments. Since this procedure is repeated after each time step, Fink 
and Soh can decrease the distances between points as desired to resolve the rollup of a 
finite vortex sheet. In addition to maintaining the spatial resolution, the method 
reduces the singular behavior induced by a vortex on neighboring points alorg the 
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sheet. However, singular behavior will appear when two parts of the vortex sheet 
approach each other closely as in the turns of a spiral. 

Baker [15] questioned Fink and Soh’s practice of amalgamating the inner region of 
the rollup into a single "core vortex". "Core dumping" is a simple approximation of 
the singularity at the edge of a finite vortex sheet. In order to analyze its effects, 
Baker used Fink and Soh’s method for describing the wake trailed by a ring wing, a 
circular vortex sheet with sinusoidal vorticity distribution. Upon observing the failure 
of Fink and Soh’s method in dealing with the double-branched rollup associated with 
the ring wing, Baker concluded that the strong core vortex must be responsible for the 
smooth rollup observed. 

Hoeijmakers and Vaatstra [16] improved Fink and Soh’s technique in several ways. 
First, they introduced the use of u sophisticated splining technique for dividing the 
aheet into panels with small curvature and polynomial vorticity distribution, 
charctcrized as a second-order panel method. Second, Pullin’s generalized similarity 
solution for rollup [ll] is invoked to make core dumping acceptable. Since 'Pullin’s 
analysis also applies to the double-branched rollup, the vortex sheet trailed by a ring 
wing can be made to rollup just as smoothly as that for a finite wing. 

Iligdon and Pozrikidis [17] introduced two-dimentional panels which are circular arcs 
with a polynomial circulation distribution. Since vortex dumping was not used, only 
closed and infinite vortex sheets were analyzed. Because the panels provide a 
continuous representation of the sheet, there is no restriction in width of the panels. 
By rising smaller panels to increase local resolution, double-branched rollup can be 
described without the use of core dumping. The existence of a panel size dependent 
instability was mentioned, but the nature of the instability was not described. 

1.4. The Motivation for this Study 

The aim of this thesis is to provide a foundation for a vortex Bhcet technique which 
offers a more accurate simulation of a thin, axisymmctric shear layer. The numerical 
techniques similar to those developed for point-vortex methods will be adopted for flat 
sections of axisymmctric vortex sheet, giving us a simple panel method for 
axisymmctric geometries. A careful, systematic study of vortex sheets represented by 
flat panels has been performed to explore the strengths and weaknesses of this 


approach. Unlike some studies of vortex sheet models, the accuracy of the mode) will 
be stressed more than the efficiency of the computational technique. It is hoped that 
much of the artificial behavior, such ns that introduced by a rigid wake analysis, can 
be removed from the simulation of vortex sheet. 

The axisymmetric panel method is obtained through the extension of the two 
dimensional panel method to axisymmetric geometries. The axisymmetric geometry is 
unique for being the only threedimcntional geometry with only two-coordinates. 
Thus, unlike the general three-dimensional flow, an axisymmetric flow can be described 
by a stream function. Fortunately, many interesting fluid dynamics phenomena exhibit 
ardsymmetry. The circular jet and the buoyant plume arc two examples of 
axisymmetric flows which can be studied easily in the laboratory. A problem of 
special engineering interest is the axisymmetric equivalent of West water’s work, the 
roll-up of the wake-sheets generated by a helicopter in hover. The roll-up of rotor 
wake is extremely important to helicopter performance because, unlike the wake trailed 
by airplane wing, the wake remains close to the rotor, significantly affecting the rotor 
performance and acoustics. 
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Chapter 2 

Mathematical Formulation of Vortex Elements 


2,1. A Review of Concepts Regarding Vorticity 

Vorticity, w, is a local property of the fluid, defined as the curl of the local 
velocity, v, 

w(x) - V x £ (2.1) 

This vector quantity describes the rotation of an infinitesmal element of fluid. A 
vortex line is a curve which is everywhere tangent to the local vorticity vector. To 
be consistent with the rotational aspect of vorticity, a vortex line must end at the 
boundaries of the fluid or form a closed curve and may never intersect itself or 
another vortex line. 

In real flows, vorticity is produced as a sheet of parallel vortex lines known as a 
vortex sheet. In invicid flow, the vortex sheet remains infinitesmaily thin and defines 
a stepwise jump in velocity tangent to the sheet The difference in velocity across the 
sheet defines the "strength” at a point on the sheet. The vortex sheet usually rolls-up 
to form what can be described as a vortex filament. In a real fluid, the vorticity 
within the vortex filament will become smoothly distributed by viscous diffusion 
within a finite time span. The structure of an evolved vortex filament is similar to 
the asymptotic limit for an invicid vortex filament with an infinite number of layers 
of vortex sheet rolled up into a region of finite cross section. In invicid analyses, a 
vortex filament with a smooth distribution of vorticity is commonly used as a model 
for a tightly wound section of a vortex sheet. 

In potential flows, the circulation 


T « ^ v . ds 


( 2 . 2 ) 


» a conserved qmmtity. This Allows the flow Associated with a complex vortex 
system to be determined by summing the contributions from a finite numljer of vortex 
elements which make up the system. The Amount of circulation in each element must 
equal the circulation In the section of the vortex system the element replaces. By 
replacing a vortex system with a set of known vortex elements with a similar 
distribution of vorticity, the flow of material points in a complex vortex system can 
be deduced. 

2.2. Two-Dimensional Vortex Elements 


2.2.1. The Point Vortex 


The two-dimensional form of the line vortex is the point vortex. The velocity 
induced by a point vortex can be derived from the conservation of circulation, (2.2), 
around a point: 


{23) 


oiknulhgi 


radial 


-r 

2itr 

0 


{2.4) 


where t is the distance from the point vortex. The velocity induced by a point 
vortex is singular when the distance from the vortex, r, equals zero. From the 
velocity expression, the stream function of a point vortex is 


-r 

tf(r) ~ — In r 
2ir 

or in Cartesian coordinates centered on the vortex. 


(2.5) 


itfx. y) ~ — In(jc 2 +y 2 ) (2-6) 

4 rr 

The two-dimensional form of a vortex filament is a vortex with finite core. If the 
cross-section of the core is circular, the flow outside the core is equivalent to the flow 
induced by a point vortex of equal strength at the core’s center. However, unlike the 
point vortex, the velocity induced by a vortex core is non-singular since the vorticity 
is distributed over a finite area. A Rankinc vortex, an example of a vortex core, 
consists of a cylinder of fluid in solid body rotation surrounded by an irrotational 


2.2.2* The Vortex Panel 


A vortex panel is ft segment of a two dimensional vortex sheet. It is equivalent in 
three-dimensions to ft strip of vortex sheet composed of parallel, straight vortex lines. 
Since invicld vortex systems would he made up of vortex sheets, accurate modelling of 
invicid vortical flow should be possible using vortex panels. 

The strength of ft vortex panel will depend on the amount and distribution of 
circulation in the panel. For this study, a panel with unit jrm distribution of 

circulation is considered. Then the gross translation of the panel can be defined at its 
mid point, its centroid of vorticity. The strength of the panel can be determined by 
its width and the ’’gradient of circulation”, K, defined mathematically as 

ar 

K — . 

dS 

The discrete form of k is obtained by the central difference scheme 


k, » 


1 H-V2 " * I-V2 
S HV2"" S \-V2 


where r\ is the circulation of the i th position on the vortex sheet and S is the arc 
length taken along the vortex sheet. Then the strength of the i th panel of width w 


is given by 


K. — 


r ^V2 " r t-Vz 
\V 


(2.7) 


Thus, a panel is equivalent to a point vortex with the same circulation "smeared” into 
the shape of a panel. 

The flow induced by the vortex panel can be obtained by representing the vortex 
lines making up the panel as point vortices of infinitesmal strengths. The velocity 
induced by the panel can then be obtained by integrating the contributions by the 
point vortices. This is can be accomplished by integrating along the panel, the effects 
of point vortices with the strength 


dr « Kds 


( 2 . 8 ) 


where k is given by (2.7). 
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In order to simplify the mathematics involved, we shall limit ourselves to flat 

vortex panels. In a two dimensional representation, the panel will be represented by a 
line segment having the Cartesian coordinates 

x(s) « x Q + s cos 6 (2.9) 

y(s) - y Q + fi slnQ (2.10) 

where the panel at (x 0 ,y 0 ) is inclined at an angle 6. 

In a previous work on panel methods Morky [18], using flat panels, carried out the 
direct integration of infinitesnial velocities to obtain the velocity field induced by a 
panel. Since the velocity integrand is singular, this involves evaluating a Cauchy 

principal value integral. However, the logarithmically singular stream function of a 
point vortex is integrable. Using liquations (2.6X2.8X2.9) and (2.10), the stream 

function is given by integrating along the vortex panel: 


\p(x,y) ln[{x — s cosd ) 2 + (y — s slnBf] ds 

4ir 

= — JU In [Q a - 2x Q 3 + s*] ds 
4tt 

« — - Q,) InQ 2 + (« + C,) lnQ t 

4lt 

a-Q a+ G 3 

+ 2Q 4 arctg — - + 2Q 4 arctg-~ 4a] 

G 9 » + y 2 

Q 4 = ycosd — xsind 
Q 3 " xcosd + yslnO 
G s * G s “ 2aQ , + 

G, = G 9 + 2aQ 3 + « 2 

where a is the half-width of the panel (a ~ Vitw). 


( 2 . 11 ) 

( 2 . 12 ) 

{2.13) 


(2.14) 

(2.15) 

(2.16) 

(2.17) 

(2.18) 


There are two methods for evaluating the velocities from the stream functions. One 
method, most suitable for Eulerian analyses with its coordinate mesh, is to take finite 
differences over the control-volumes. 'Ilxc alternate method is to evaluate the velocities, 
which arc derivatives of the stream functions, explicitly at the points of interest. 





The velocity field produced by the equations are shown in Figure 2-1. In these 
formulas, the arc-tangent, being multiple valued on the panel, models the jump in 
tangential velocity as the panel is crossed. Mathematically, the panel is a branch cut 
in the coordinate space connecting the two logarithmic singularities at the edges of the 
panels. 


2.3. Extension to Axisymmetric Geometries 


2.3.1. The Axisymmetric Problem 

An axisymmetric vortex sheet is composed of coaxial, circular vortex lines. Thus, a 
simple-minded extension tf the Rosenbead approach to sheet discretization is to replace 
point vortices by circular line vortices. However this straight-forward approach is 
unworkable because a curved line vortex is influenced by its own velocity singularity. 
As a result, a circular line vortex will translate through the fluid at infinite velocities 

[19}. One way of avoiding the infinite translation velocity involves using circular 

vortex filaments, or vortex rings, in place of the singular line vortices. This 

introduces a finite vorticity distribution which induces a finite translation velocity. 
Another approach is to use polygons made of straight segerants of line vortices in place 
of the vortex rings. A network formed by vortex polygons is known as a vortex 
lattice. Unlike the strictly axisymmetric models, the vortex polygons should be able to 
display nonaxt. symmetric deformations of the vortex sheet. Azimuthal deformations arc 
known to hasten the decay of vortex rings [20], and should have an important 
consequences for the behavior of three-dimensional vortex sheets. 

However, both of the above models of axisymmetric vortex sheet require the selection 

of a length parameter: the radius of the vortex core or the length of a polygon 
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element. Thun before any of these models cun be used accurately, it is necessary to 
gain a better understanding of the behavior of vortex sheet these models are designed 
to reproduce. 

The axisymmetric equivalent of the vortex panel is a "vortex band”, a aection of an 
axisymmetric vortex sheet. It represents a finite piece of the vortex sheet which, in 
some models of the axisymmetric vortex sheet, is replaced by vortex rings and vortex 
polygons. Although a band docs not poses* any volume, by defining the band’s 
position at the mid point of the band as was done for the tvm-dinicntional panels, the 
self-induced velocity of the band is shown to be finite. 


The stream function for circular line vortices, expressed using complete 
integrals E l (X) and I^CX), is given by Umb [21], 

«x.y) = •— (r, + r 2 ) { jF,(X) - £,(X) ) 

2ir 


r, = Least distance between (x,y) and the vortex 
r 2 s Greatest distance between (x,y) and the vortex 


elliptic 


(2.23) 


(2.24) 


To obtain the flow for ft vortex panel, this stream function must be integrated across 
"he panel, resulting in a difficult double-integration of singular elliptic functions. 


2.3.2. The Method of Matched Asymptotic Expansion 

If the width of the vortex ban,’ is small relative to its radius, an asymptotic 
expansion for the stream function can be obtained. The approach used here for the 
vortex panel, known as the method of matched asymptotic expansion, was also used In 
a study of curved vortex filaments [22], Ibis method derives an uniformly valid 
solution which asymptotically approaches known behavior of the actual stream 
function. Thus, it is important to hcuiistically deduce what is expected from the 
stream function. We can expect the character of the stream function away from the 
hnnri tn (innfuarh that nf the stream function for a circular line vortex. Thus, the 

•****” -jti 

circular line vortex stream function will be used as the outer-solution. Near the band, 
the stream function should resemble the two-dimensional solution for a vortex panel. 
However, the inner-solution must include effects due to curvature in the axisymmetric 
geometry. 


The coordinate systems used in the following analysis are shown in figure 2 2, 

2.3.2. 1. The Outer Solution 

The velocity away from the band is derived from (2.23) which gives 

"... - j ly < + ’■) i - *,w > > <m 

' 4 ^ 

v «* m ~~ (>\ + »'*) / " #,<*) 7 j (2.26) 

where H is the radius of the vortex band. 

2.3.2.2. The Inner-Solution 

The axisynunetric velocity field, for a given axisymmctric stream function \p, is 
defined to be 


u[r,z) « 


zi 

r dz 


v(r.*) ~ 


1 d\l/ 
r 


The irrotationality of the potential flow gives 


(2.27) 

(2.26) 


du dv 

— ~ B 0 

dz dr 

Which is the same as that for twodimcnsional flow, 
(2.27X2.28) and (2.29) gives the governing equation 


(2.29) 

w * 0. Combining Equations 
for the axisynunetric stream 


functions: 


d*+ ^ d 7 j/ 1 d* __ 

dr 2 dz 2 r dr 

This equation is the axbymmetric equivalent 
equation 


( 2 . 30 ) 

of the two-dimensional stream function 


d 2 $ 

dx 2 


+ —2 " s* o. 


dy 2 


(231) 


It is now necessary to switch to a panefcentcrcd coordinate system (x,y), originating at 
the center of the vortex band of radius R. 
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OH + OH ^ i (fy 

Ox 2 Oy 2 ~ Mix fc m a (2,32) 


Non tllmpjiNlonnllzIng with resect to the outer length scale, radius JR, gives 

gy f OH 1 0$ 

Of l+x Ox " °‘ (*M) 

coordinates can now be rescaled by a small parameter, the curvature, which is 
* * » 'where a is the halt-width of the band. Rescaling (2,33) by e gives 

OH + rJ H e_ 

r)y* J+e* a* * °‘ (Wf) 

with the understanding that (x,y) are non- dimensional, stretched coordinates. 

Tiie expanded inner stream function. 


* ' * + **■ * '**. • • • 0.35) 

” "* d “ tmCSS eff ““ of '• *» to the perturbation of the two- 

dtmcnsional atream function. *, In the .hove Krira , ttrms w ,» M hlve ^ 

included with torn* <X.-» ine). Ualng thi. eapanaion, the governing et,u«ion0.34) 
can be broken up into equations for each order in the expansion: 


(236) 

(237) 


’ — “ jjuvcuung equation tor the two dimensional stream 

funetton(2.31). The autaquent equation, for the highor order tern, of the inner .tream 
unction are all Puisrena Equations with the non-homogeneoua term, involving the 

aolutiona of lowemrder equation, Thu, repetitive evaluation uf the Maaun’a Equation 
can be used to obtain the desired degree of accuracy. 

By aasumlng the effect nf curvature to he entail, we will only need to evaluate the 

tat term to the eapansion. The particular solution for the firatmler I'om.n’s 
liquation is 


0(1) • 

OH* + 0H e 

- 0 


ox 2 

dy 2 

<Xe)* 

0 2 \p i + 0H X 



ax 2 

dy 2 

Ox 




x W 


( 238 ) 


Then, using (2,13) for the zeroth order term of the expansion, the inner stream 
function is 


K m (> + - 2 *x) + V'hwtvi' (2. 99) 

where 0 horno is some stream function satisfying the homogeneous equation, which is 
determined by matching the inner end outer solutions as shown in the following 
section. 

Using notations defined by (2. 14X2. 15X2.16X2. 17) and (2.18), The inner solution for 
the velocity field is 

“in " (1 + , c x) [sin $ Iny?- + 2cos Q (orctg--™ + arctg -—)] (2.40) 

4n A «, w* 

k ; Q.. «+Q 

v m “ U + 2 ex)l-cos0ln ( yi ■2slnd(arctg^^4arctg”-rj^)J (2.41) 

+ \ < /(« “ Q 3 ) lnQ z + (a + e 3 ) lnQ y - 4a 
+ 2Q a (arctg — ^ + arctg—^ . 

2.3.2.3, The Inner/Outer Matching Solution 

The vortex band inner stream function solution (2.39) must be matched to the outer- 
solution, which is approximated by the circular line vortex solution (2.23). First, the 
circulations of the perturbed vortex panel and the circular line must be equal. Then 
it is necessary to equate the outer-limit of the inner solution to the inner-limit of the 
outer-solution using a matching criteria. The "matching solution’*, the intermediate 
limit solution obtained using the matching criteria, will subsquently be subtracted from 
the sum of the inner and outer solutions to produce the uniformly valid asymptotic 
solution. A discussion of the method of matched asymptotic analysis can be found in 
Reference [23], 

The inner-limit of the outer solution is found by expanding the elliptic integrals in 
(2.23) for small complementary modulus, 1 - X 8 (2.24), as shown by Jahnkc and 
Fmde [24] 


dr' 

*oul 



2 7T 


+ 


X 

2 ft 


+ ln 8R ^ + + <J ( f 2ln *) 


This result was also obtained by Tung and Ting [25], 


(2.42) 


The outer limit ol' the inner solution is found at the limit, r ■,*>«, where the panel 
width and inclination lwcome unimportant. Therefore, using the approximations 


C>, ^ x v + f * r* 

Q r- x ? + f K 


hi 


0 


«C> 3 „ 

arctff -zr • + on% , 4 *=- ; «rctg #1 — ^ 0 

**4 


e. 






wc can get the outer-limit of the inner-solution, 


h° 

'In 


-Ka r. x . 

7 ( ! + 5*> 


neglecting the higher order terms. 


(2.43) 


Since the panels span the width 2a, the definition for k(2.7) 


r p ft n e l " 2 * * 


(2.44) 


is substituted into (2.42). The complete intermediate stream function for a vortex band 
is produced the matching criteria, ^ horoo . which equates the two intermediate limits 
(2.42) and (2.43). Neglecting the higher-order terms, the matching criteria is added to 
the inner-limit of the outer solution (2.42) 


— Ka , , x 

K - — < 2 + 

IT 


2R * ^SR ^ + 2R^ * ^ hu ™‘ 


(2.45) 


The matching criteria, ^ ho mo 48 obtilincd b y equating the inner-limit of the outer 
solution (2.45) to the outer-limit of the inner solution (2.43) 


K . { ,n~ (J + 

U * §R + ln M V + S>' + *- 


i r 

’Ka 


Then the matching criteria ia found to be 


~ f* + £ + ln l < J + £»■ 


(2.46) 


(2.47) 


'horo _ 7R 8 ' 2,R 

and the intermediate stream function is established, as shown by liquation (2.46), to be 
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-KH . X . r y 

w ' lt * 2*> '"«■ 


(2.4H) 


2. 3.2.4. The Composite, Uii|t(i>’inly V(ilM Solution 

The complete, uniformly- valid stream function for the axlsymmetrlc vortex Irnml In a 
composite of the Inner solution, the outer solution, \ p^ n , and the Intermediate 

(matching) solution, \p^ n . The complete solution can l« formulated to Iks 

Mx.y) « * )rt + \p m - {2.49) 

where, 

^in • -Inner Solution (2,? 9) 

= * O + 2*> K 

\p i0 r: 7’wo -Dimensional Panel Solution (2,13) 

“—(«*- Q ) lnQ 2 + (a + Q 3 ) InQ 
4it 

a-Cs «+fi, 

+ 2Q i arctg-- + 2Q li orctg— - 4a 

^ sut ~ Outer Solution (2.23) 

* — (r,+r 2 ) { jP,(X) - #,(X) ) 
it 

^ |o ™ Inner/Outer Matching Solution (2.45). 

"™Ka , T \ . x , , 

« ~ ii ^ In — ( 1 + -r r) J 

It a 2K 

using notations introduced earlier. Notice that ^ lu and have been 

redimensiomlizcd with respect to the outer-dimension, R. 


The above stream function is antisymmetric, thus its velocity field is defined to be 


«(x,y) « 

v(x.y) - 


- 1 

l)y 

1 dty 
* dx 


However, the inner-solution, \p in , and the intermediate-solution, arc valid only for 
the inner region, x « R. Thus, the inner -region velocity field, using Rrx r- R, becomes 
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( j 

i y <*» - 

~ v “ “ s t ( *“ - *J- 


while the outer velocity field remains 


j \ * 

u( X ,y) * — ~~ 

v(jc.y) 


7 ^ 
Sy 


ii+x 

Thus, the velocity field is given by a composite solution of the form 


u(x,y) = « |n + « ou( - « to 
v(x.y) - v (n + v cu , - v h 


(2J0) 

(2Si) 


where, 


u * 

In 


4it 

K 


„ Q 2 «“V, an-y. 

2 ^) l~sln 6 + 2cos d (arctg~~y— + arctg 

*1 ^4 W4 

( J + 2j^ ^"COf<9in™-2jlne(arcf^^ + arcfg^™-); 

™ /(a~£> 3 ) i«G 2 + (a+G 3 ) inQ, - <te 

a ”G 3 o+Q, 

+ 2Q 4 arctg + 2G 4 «rcf g~~r~] 

«4 y« 

Hy + yKR-li) + (r,+r,)?5^-5; 

r ’ r ’ 2 dy 


«“G, 


ou( 


K a 

rtf/Hx) 
•na 


I '8 
* . *+21? , 


*{£+*) r i r * dy 


I '8 

* , y 


, * * y 


V MS . 
to 


Ka 


x , X 






■2 2Jk"'8R J ‘ 


The relationship between part? of the complete solution is shown by Figure 2-3, 
showing the contributions to the velocity profile taken radially from the cylindrical 
coordinate axis, denoted by the center-line, through a vortex band « 0). The 
uniformly -valid, composite solution for the velocity induced by the vortex panel, 
v comp’ ^ denoted by crosses. I'rominent features in the plot arc, 


23 


1) the 1/r type singularity of the outer-aolution at the center of the 
panel Cx « 0, x x « 0) matched by the intermediate-solution, 

2) the logarithmic singularities of the inner-solution at the panel edges. 

The composite solution retains the outer solution's 0(e Ine) portion on the axis of 
symmetry while smoothly blending in the two-dimensional inner-solution in place of 
the singularity at r^ -0 before decaying off faster than the two-dimensional solution. 
Notice that the inner-solution is displaced, representing the speed of the self-induced 
translation of the vortex panel. 


2.3.3. The Self-Induced Velocity 


An important result of the preceding analysis is the self-induced velocity of a vortex 
band. The motion of the vortex band is defined at the mid-point of the band. Thus, 
using Equation (2.51), the axial speed of the panel center is found to be 


bond 


* — Mx-O.y-O) 
dx 

Ka 8R J 

= — {in— + -J. 
2irR a 2 


(2.52) 


The radial velocity at the panel center is zero for any orientation. Thus, for the 
present order of accuracy, 0(e), the self -induced translation of a vortex band, defined 
at the center of the band, is purely axial and independent of the panel orientation. 


Ry applying (2.44) to equate the circulations, this mult can be compared to the self- 
induced velocity of a uniform vorticity vortex ring, where ™ is constan;, [21, 22] 


ring 


Ka HR 1 

{ In — - 

2nR r 0 4 

Thus, the core radius of an "uniform* 


vortex ting must be 


(2.53) 


R s * {2M) 

0.4724 ~ 

A\ 

in order to translate at the same velocity as a vortex panel of the same circulation 
and radius. Since the self induced velocity reflects the energy of the vortex element, 
the ratio pursued by equating the known energy of a vortex disk, with the energy 
numerically integrated for a set of equally spaced vortex rings. The core diameter of 


the vortex rings giving the same amount of energy is obtained numerically and the 
ratio is the quotient of this value over the ring spacing. By increasing the number of 
vortex rings, the ratio was found to converge to about 0,50 [261. 
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figure 2-1: Velocity Vectors Associated with a Vortex Panel 
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Chapter 3 

Numerical Behavior of Vortex Elements 

The FORTRAN subroutines for calculating the velocities induced by different vortex 
elements are presented in Appendix A. In this section, the differences between vortex 
elements will be explored. 

3.1. The Velocity Profiles of Two-Dimensional Elements 

Figures 3-1 and 3-2 show velocity profiles for vortex panels. The velocities are 
perpendicular to the panel -centered coordinate axes because of the panel’s uniform 
circulation distribution. The most important feature is the lack of a strong sigularity 
in the panel velocity profile. However, the normal velocity has a logarithmic 
sigularity at the edges of the panel and no singularity in the tangential velocity, the 
latter being double-valued but finite on the panel itself. Away from the panel, the 
velocity profile quickly approaches the point vortex profile. The difference in the 
effect on the adjacent vortex elements is not obvious in this representation. 

3.2. The Velocity Profiles of Axisymmetric Elements 

Figures 3-3, 3-4, 3-5, and 3-6 show the velocity profiles along the panel axes of vortex 
bands with inclinations of 0 and 90 degrees. The centerline to the right denotes the 
location of the axis of symmetry with respect to the width of the panels. The self- 
induced translation velocity of the vortex band appears as the non-zero axial velocity 
at the midpoint of the panel. 
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3.3. Integration of Vortices In the Form of a Vortex Panel 

An a simple test of integrating vortex elements to simulate ft vortex sheet, vortex 
elements were arranged to form ft vortex panel. The flow produced by ft vortex panel 
of uniform circulation is described completely by Equation (2.13) and provides n 
baseline for the flow induced by the integration of vortex elements. Figure 3-7 Bhows 
the velocity vectors of flow around the baseline panel. The figures 3 8 and 3 9 
shows the flow induced by four and ten point vortices, respectively, arranged to 
simulate a vortex panel. The accuracy in reproducing the baseline flow demonstrates 
the gain in accuracy achieved by integrating more point vortices. Of special 
importance is the velocity induced on the panel represented by the point vortices. The 
velocities induced on the point vortices are shown along the dashed line representing 
the vortex panel the point vortices replaces. For four point vortices, the singular 

nature dominates the velocity near the panel. For ten point vortices, the velocity field 
has a closer resemblence to the baseline flow. In comparison, Figures 3-10 and 3-11 
show the flow field induced by two and ten vortex panels respectively. Unlike the 
previous plots, the velocity along the integrated panel is determined along the 
coordinate grid, they do not necessarily represent the velocities at the centers of the 
panels. The uniform distribution of velocity along the integrated panel shows the 
effective cancellation of the logarithmic singular behavior at the panel edges. The 

overall velocity fields are identical, indicating the validity of replacing a vortex sheet 
with an integration of vortex panels. 
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Figure 3-7: Baseline Flow: a single vortex panel spanning x - ±1 
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Figure 3 8: Point Vortex Flaw: 4 vortices 
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Figxire 3-9: Point Vortex How: 10 vortices 
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Flgnrc 3-10: Vortex Panel Plow: 2 panels 
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Figure 3-11: Vortex Panel Plow: 10 panels 



Chapter 4 

Discrete Numerical Approximation of a Vortex Sheet 


4.1. Th«> Uniform Downwash Solution for Dlaerotized Vortex Shoots 

The simulation of tiic vortex sheet trailed by a wing with an elliptic loading has 
been the most popular application for numerical vortex sheet methods. This popularity 
is not only due to the practical significance of the elliptic loading but also because the 
loading is predicted to generate only one pair of roll-up spirals. 


An elliptical ly loaded lifting line, with semi-span b, poetesses the circulation 
distribution of the form 


n*) = r, / 1 - ~ f». [ 4 . 1 ) 

A straight lifting line with this distribution will generate a uniform downwash and a 
flat vortex sheet into wliich the bound vorticity is lost as one proceeds outward from 
mid span. Thus the strength of the sheet is given by 


k(jc) 


dr 


dx 

r . ? 11 


v2 


6 2 


1 V* 


( 4 . 2 ) 


Because of the singularities at the tips, the sheet will immediately roll-up into 
exponential spirals. The rolling-up of the sheet relieves the discontinuities in the 
vertical velocity profile. 


Numerical simulations of the wake sheet dynamics usually ignore the fluid dynamics 
of the vortex sheet formation. For purposes of numerical analysis, it is important to 
know the initial geometry of the sheet and the flow associated with it. In many 
Lagrangian approaches to roll up, an initial distribution of discrete vortices is selected 


without Checking its ao , ;u r ncy In modelling the initial state of the vortex sheet, and 
tlic merits of the analysis are judged by the quality of the resulting roll up. 

A flat vortex sheet which indtuxs a velocity field which is constant along the sheet 
and singular at the tip is not adequate as the initial condition for a roll up 
calculation. For a realistic model, the flat vortex sheet model can only be relied for 
an "outer-solution'’, valid away from the singular discontinuity at the edge of the 
sheet Then an "inner solution” at the edge in needed to give the effect of the sheet 
roll-up which is expected to fori*'., immediately upon the formation of the vortex sheet. 
Thus, a realistic simulation of the wing wake must begin with some sort of rcll-np to 
replace the tip singularity. The same conclusion appears to have been reached 
empirically by lloeijmakeis and Vsastra [16] 

4.1.1. Point Vortex Representations of the Vortex Sheet 

For Lagrangian simulation of vortex sheets, the initial downwash profile is 
represented by the velocities induced on each vortex element. Surprisingly, past 
investigators did not evaluate or comment upon the accuracy of the vortex sheet molel 
to duplicate the initial condition. The only published result which gives the initial 
downwash profile is West water’s pioneering work. In Westwater’s initial vortieity 

distribution, an error in the position of three point vortices produced a "kink" in the 
initial downwash profile as shown by Figure 4-2. Tire total circulation contained by 
the sheet is 1.0, resulting in the uniform downwash result of 0.5 for a continuous, 
flat wake vortex sheet. Although the kink is apparent at the beginning of the 
tabulated results [12], Westwater proceeded with the tedious manual calculations 
without correcting the error. Unfortunately, even though the initial downwash profile 
can indicate the accuracy of the vortex sheet model and the validity of the subsequent 
foll-up, such results have not been published with the results of other roll-up 
calculations. 

There are two approaches to representing a vortex sheet using point vortices. 
Westwater introduced the method of dividing the sheet into strips of equal circulation 
and replacing them by a point vortex at the centroid of vorticity of the strip. The 
other method, used successfully by Fink and Soh, replaces strips of equal width with 
a point vortex of equivalent circulation at the center of the strip. Westward's method 


* * 
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WM *'° imlar for * 101 ‘2 timB " iuc « ^ concentrates the vortices. «iv! *1 ko served ns tracer 
particles, where mat action i* expected. Pink am! Soh, on tire other hand, designed 
their method to approximate an integration of flK panel* The downwash profile 
calculated using these two methods are shown in Inures 4 3 and 4 4. Both methods 

share the smoot «i, yet large, deviation from the uniform down wash value of 0.3, 

which is improved by increasing the number of vortices used. The exaggerated' 

smoothing of the velocity discontinuity appears to be responsible for the extremely 
amooth roll-up Obtained using point vortices. Figure 4 4 displays a secondary deviation 
by the velocity of the inner-must point vortex of the distribution. This error is 
caused by the proximity to the mirror image of the vortex which, having opposite 
circulation, adds a significant contributes on the downward velocity of the inner-most 
point vortex. In Figure 4-3, where there Is more separation between the inner-most 
vortex and its adjacent vortices, thus no deviation is observed. Overall, the deviations 
from the predicted value of the uniform downwash are surprisingly large. Point 
vortex correction derived by Van der Vooren 1 , as presented by Moore [4 is not 

applicable for the equal space discretisation. Application of Van dcr Vooren’s correction 
to the equal strength discretization did not improve the calculated downwash profile. 


Instead of using the velocities induced on the point vortices to calculate the motion 
of the vortex sheet, the downwash can be defined to be the velocities obtained at 
points half way between equally spaced vortices. This is the two-dimensional 


equivalent of the quadrilateral vortex-lattice (or ’’dipole”) method which is often used 
for three-dimensional flow simulations. Figure 4-5 shows the velocities at points half- 
way between point vortices representing an elliptkally loaded wing wake. The 
vertices are equally spaced as in Figure 4-4. The downwash profile obtained by this 
method differs substantially from the uniform lifting line downwash (0.5) and displays 
a character different from Figures 4-3 and 4-4. The "dip” in the downwash at the 

end of the profile is the result of the large step increase in the vortex strength at the 
tip. 


./van dor Vooren'a correction involvi 
distribution, resulting in an addition 
to the velocity induced on each vortex. 


the gradient of the vorticity 
point vortex-like contribution 
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4.11.2. Vortex Panel Representation of the Vortex Sheet 

The downwash profiler corresponding to Figures 4 3 end 44 for integrations of 
vortex ptjacl" are shown in Figures 4-6 and 4-7 respectively. The figures show that 
an integration of panels will produce much better approximation of the downwash 
profile expected for a flat vortex sheet than possible using the some number of point 
vottiee.**. Near the tip, the deviations from the uniform downwash result reflects the 
large step increase in the discretized distribution of the vortex sheet strength. The 
equal circulation panel disttibution, Figure 4-6, produces less panels where the sheet i« 
weak, i exulting in the inaccurate velocity on the left-most panel. Since the velocity 
of a vortex panel is obtained at tlx center of the panel, which are half-way between 
the logarithmic singularities at the panel edges, features found in the velocities obtained 
between point vortices axe expected. The only observed common feature is the 'kink” 
in the profile nt the next-tcKhe-lust panel also found in Figure 4-7. However, the 
kink in the velocity induced by panels is less revere than that for the velocities 
midway between point vortices. This reflects the difference in degree of singularity 
possessed by the two vortex elements. 

The downwush profile for the integration of equal width vortex panels closely 
resembles the uniform downwash produced by an elliptically loaded lifting line except 
at the tip. The discontinuous nature of the downwash inhibits the smooth roll-up of 
a sheet modelled using only vortex panels. If the downwash induced on the panels 
(disregarding the last panel at the tip) is taken to represent the outer-solution, then an 
inner-solution must be superposed to represent the roll up at the tip. The tip roLl-up 
can be represented by replacing the lest panel on the sheet with a vortex core of the 
same strength, positioned in accordance to the Kaden’s similarity solution (See Appendix 
B for details). Figure 4-8 shows the resulting smooth downwash profile. The right- 
most point in the profile is the downwash on the last panel and its "relief” in the 
downwash profile is significant; the velocity of the tip roll-up core is positive and lies 
outside of the plot. 

The downwash profile for the integiation of equal strength vortex panels appears to 
show the effect of a point vortex near the tip, the virtual core- The effect of such a 
core can be demonstrated plotting the downwash profile calculated with a point vortex 
of the opposite sign at the tip to cancel out the virtual cor* t st (Figure 4-9). The 


shape of the downwash profile near the tip i* not sensitive to the number cf panels, 
suggesting some sort of similarity solution at the tip. The velocities induced on the 

last three panels by the integration of 5 to 100 panels were used to determine the 

strength snd position of the point vortex needed to cancel out the deviation from the 
uniform downwash value. The results, plotted in Figures 4 10 and 4-11, shown that 
the vortex must be positioned 0.204 times the half-width of the last panel from the 
edge of the last panel and contain about 0.0547 times the circulation contained in one 
panel. 

4.1.3. Vortex Band Representation of the Vortex Sheet 

The axisymmetric equivalent of the wake vortex sheet of an clliptically loaded 
lifting line is the disk of axisymmetric vortex sheet w. l represents the translation of 
a disk in a fluid at rest. The first theoretical study of this flow was made by GJ. 
Taylor [271. The vorticity distribution necessary to obtain a uniform velocity profile 
on the disk of vorticity is equivalent to liquation (4.2). The only difference is the 

factor of ^ greater downwash induced by the same amount of total circulation. This 

means that using the same vorticity didtribution as before, we expect the uniform 
velocity to have the value of ~ - 0.7854. 

The straightforward discretization of the axisymmetric vortex sheet into circular line 
vortices will not be meaningful due to the infinite self-induced velocities of the 
vortices. However, if the circular line vortices are replaced by vortex rings with a 
finite core dimension r Q , then the self-induced velocity of the vortices will be 0( olno\ 
where a is the ratio between the core dimension and the ring’s radius. In addition, 
since the core represents a section of the sheet of approximately the same size, the seif - 
induced velocity is expected to be ita R lno) when compared to the 0(1) velocity 
induced by the total disk of vorticity. Thus, If O Is small, the effect of the self- 
induced velocity is expected to become negligible. If we hypothetically allow the 

infinite self-induced velocities of the circular line vortices to be neglected, then the 
discretized disk of vorticity will induce the downwash profile shown by Figure 4-12. 
The downwash profile for a set equally spaced circular line vortices fc surprisingly 
uniform and near the correct value, thus demonstrating the insignificance of the self- 
induced velocity. 
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If the moving disk is dincratlied by vortex bonds, the seif-induced velocities of the 
bends may be included. The resulting downwosh profile is shown in Figure 4 13, 
liiccpt for the magnitude of th* velocities, the profile displays the same features found 
in the wing wake profile, However, it should be noted that the accuracy of the 
velocities docs not uniformly improve with the increase in the number of vortex 
bends. Figure *14 shown the use of a tip-core to approximate the inner-solution 
representing the roll-up of the tip (Sec Appendix B for details). The shape of the 
resulting down wash profile near tire tip does not agree with the two-dimensional 
equivalent produced by Kadcn’s similarity solution (Figure 4-8). This m:.y be an 
indication ths.t the core is positioned too far away from the rest of the sheet in the 
approximation used lo model the tip region of an axisymmetric vortex sheet. Since the 
approximation is not a similarity solution for the geometry of the tip roll-up, the 
geometrical error will be reflected in the distance between the edge of the last panel 
and the tip-core. 

4.2. Summuary 

Tire calculated velocities of the vortex elements can be used to test the accuracy of 
the vortex sheet model if the velocities of the actual vortex sheet are known. For 
elliptical loading cases discussed above, the velocities for each section of the vortex 
sheet can be predicted. These cases also represent the initial conditions for interesting 
problems in vortex sheet dynamics. The above results show the accuracy of equal 
width panels in reproducing the flow field away from the singular tip of the vortex 
sheet. A near-field solution for the fiow around the tip of the vortex sheet can then 
be approximated using a core of vorticity. The importance of the tip core is also 
suggested by the virtual core effect in the downwash calculated for the equal 
circulation segmentation of the same vortex sheet. The incorporation of a tip core in 
the equal width panel method represents a ctraightforward modelling of the vortex 
sheet with results which arc more realistic than the results obtained using point 
vortices alone. 
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Figure 4-1; Vertical velocity profile for elliptic loading 
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Figure 4-3: Downwash Profile: Integration of equal strength vortices 



Figure 4-4: 


Downwash Profile: Integration of equally spaced vortices 
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Figure 4-7: Downwash Profile: Integration of equal width vortex panels 
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Figure 4-8: Downwash Profile: Kqual width panels with roll up model 
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Downwash Profile: The effect of virtual core 
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Figure 4-10: 


Virtual Core location: Normalized w.r.t. the tip panel 
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Figure 4-11: Virtual-Core Strength: Normalized w.r.t. the tip panel 
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Figure 442: Moving Disk: Circular line vortex representation 








Chapter 5 

Numerical Simulation of Vortex Sheet Dynamics 

5.1. The Numerical Roll-Up of Vortex Sheets 

The roll-up of two-dimensional and axisymmetric vortex sheets were simulated using 
vortex panels. The changing shape of the vortex sheet is determined in a Lagrangiun 
fashion by following the vortex panels throughout cycles of numerical, time integration. 
Accuracy of the method is maintained by reconstructing and rcdiscretizing the sheet 
time integration cycle in the manner prescribed by Fink and Sob [131 A first-order 
Runge-Kutta scheme is used for each time integration cycle. The main difference 
between the panel method and the point vortex method is the amount of data which 
must be handled to describe the sheet. Since the orientation and the width for each 
panel must be recorded, the amount of geometrical data is effectively doubled. The 
numerical scheme for following the dynamics of the vortex sheet is summarized by 
the flow chart in Figure 5-1. 

5.1.1. The Time-Integration Scheme 

The centers of the vortex panels arc translated using the Runge-Kutta scheme during 
each cycle of time integration. A vortical flow system is characterized by gradients in 
velocity resulting in curved trajectories of the converted vortices. The Runge-Kutta 
scheme is generally more efficient in describing curved trajectories than the simpler 
Euler scheme [28J, Euler scheme translates the control points for an short increment 
in time along a straight trajectory based on the velocity at the beginning of each 
cycle. Runge-Kutta scheme, on the other hand, uses more than one set of velocity 
data to approximate a curved trajectory during each time integration cycle. If only 
one additional velocity input is used, the Runge-Kutta scheme is Tefered to as being of 
the first order. The panel method, unlike the point vortex methods, requires the size 
and the orientation for each panel in order to calculate meaningful velocities. Since 
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only the cnitcir of the panels me tracked by the time integration, the entire sheet 

must Ik; xrrniiNtructed in order to obtain the «lze and the orientation of the panels. 

However, n rigorous reeonsti ruction of the sheet is highly computation Intensive. 
Prohnbly for this reason, Hoeljmuktw anti Vantatra |U>] used the Euler scheme with 
their panel method, while l ink and Muh (14] were able to use a Kunge Kutta scheme 
with their improved point vortex scheme. 

in the fits* order Kunge Kuua scheme used in tills study, u simple approximate 

reconstruction of the lute* mediate vortex sheet geometry is used to provide additional 
velocity data. At the beginning of tire Integration cycle, the vortex panels are 

arranged by the rcdiscretization scheme to form a continuous vortex sheet and the 
initial set of velocities are obtained. Then the intermediate positions of the panel 
centers are derived using the Euler scheme. To obtain, the intermediate velocity, the 
changes in the orientaion and size of the panel due to the Initial deformation of the 
sheet must be derived. The intermediate panels arc oriented parallel to the line 

spanning the positions of the two adjacent panels center's (figure 5~2). The width of 
the intermediate panel is taken to be half the distance between the adjacent, panel 
centers. The new positions of the panel centers are obtained by averaging the initial 
and intermediate velocities and by translating the panels by the correc'.ed velocity for 
twice the time increment used to obtain the intermediate solution. 


5.1.2. The Reconstruction and Rcdiscretization of the Vortex Sheet 

The reconstruction and rediscretization of the vortex sheet are performed after each 
cycle of time integration. The rediscretization procedure allows the vortex sheet to be 
reproduced by vortex panels of equal widths. This scheme also helps to maintain the 
accuracy of the panel representation of the vortex sheet and helps dampen short 
wavelength instabilities. The damping of the vortex sheet is achieved by reconstructed 
the vortex sheet using only the positions of the panel centers. A parabolic spline, 
selected in ovder to rule out non-physical wiggles between the control points, is fitted 
to represent the vortex sheet (figure S 3). At the same time, a linear interpolation is 
used to represent the distribution of circulation in the sheet. The new vortex panels 
having equal widths, are made to span segments of the vortex sheet (figure 5-4). The 
circulation contained in the segments of the vortex sheet is transfered to the panels 
spanning them. Although this procedure reduces the effective radius of curvature of 


the vortex sheet by shifting the panels inside a curved spline representing the vortex 
sheet, the effect can 1* minimized by decreasing the size of the panels. In sny case, 
because straight segments are being used to represent a curve, email panel width is 
necessary to maintain the geometrical accuracy of the panel representation. The size of 
the segments are also adjusted to maintain a certain relative distance between the 
outer most panel and the vortex representing the Inner region of the roll-up spiral. 

5.1.3. The Treatment of tho Inner Roll-Up Region 

As shown previously, the vortex panels alone can at most reproduce the outer 
solution for the flow induced by a vortex sheet. This is inadequate as an initial 
condition without a model to represent the initially infinitesmal roll-up at the end of 
the sheet. Thus, the simulation of vortex sheet rollup begins by replacing the 
singularity at the edge of the sheet with a core representing the initial roll-up of the 
sheet edge (Figuro 5 5). The placement of the tip roll-up core for the two-dimensional 
wing wake follows the result presented by von Karman of Kaden’s analysis for roll- 
up by a parabolically loaded vortex sheet. l J or the roll-up of oxisyrametric vortex 
sheets, due to the lack of a similarity solution, the a much simpler tip core is used. 

The details of the initial conditions are described in Appendix B. Such models for the 

initial condition of the vortex sheet is necessary for a realistic numerical approximation 

of the flow when using vortex panels to model the vortex sheet. 

The core is a model for a tightly rolled-up section of vortex sheet, the core center 
representing the center of the spiral. By knowing the position of the center of rollup, 
the inner-most resolvable section of the vortex sheet in the spiral can be handled more 

accurately. The spline representing the vortex sheet is derived from the panel centers, 

the location of the edge of the last panel, corresponding to the end of the sheet and 
the inner-most resolvable section of the spiral, is uncertain. (See Figure 5-6.) Since 
the roll-up should be in the form of an exponential spiral with the core at its center, 

the orientation of the last panel may be used to blend the resolved section of the 

spiral with the core. It also possible to truncate inner-most section of the resolvable 
vortex sheet and dump the circulation, momentum, and energy of the truncated section 
into the core. (Although both methods may be used concurrently, core dumping was 
not used in any of the results presented below.) Core dumping is attractive for 
numerical reasons since the size of the panels can be kept relatively small even as the 
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,hcct stretches without toeing limited by computer memory size. As vorticity is 
dumped into the core, the conservation of circulation, momentum and kinetic energy 
determines the strength, position, and size of the core respectively. 

5.2. Result's of the Roll-Up Calculations 

The numerical scheme described above has been applied to simulate the roll-up of 
two-dimensional and a asymmetric vortex sheets. All roll-up simulations were initiated 
with a model representing a small amount of roll-up already in place at the edge of 
the vortex sheet. 

5.2.1. Two-Dimensional Roll-Up 

The roll-up of a two-dimensional vortex sheet trailed by a flat, elliptically loaded 
wing is the classical test case for numerical vortex sheet methods. In order to 
simulate the geometry of the sheet, the Trefftt plane method, originally introduced to 
roll-up studies by Westwater [121 is used. Each cross-section of the wake vortex sheet 
is treated as an unsteady two-dimensional vortex sheet. Thus, as described above, the 
sheet is replaced by a series of equal width flat panel of constant vorticity. Unlike 
the axisymmetric flow system, this two-dimensional roll-up requires the program to 
account for the mirror-image contribution across the plane of symmetry. The time 
integration cycle is begun with a flat sheet at y - 0 and a tip core position^ at the 
center of Kaden’s spiral for roll-up of the outer 5% of the vortex sheet. There are 76 
panels describing the sheet initially, the number increasing as the sheet stretches; the 
size of the panel is kept constant. In this simulation, roll-up of the sheet has been 
allowed to take place without core dumping to demonstrate the robustness of the 
model. The panel widths have been kept constant. 

Figure 5-7 ehows the geometry of the sheet after every ten time integration cylw, 
each cycle representing a nondimcnsional elapsed time 2 of the the of 0.005 . The 
initial location of the roll-up core is given by a small circle, partly hidden by the 
inner regions of the spirals, at about x « 0.97 . Although the inner region of the 
spiral is obscured by the overlapping, smooth rollup of the outer region can be seen. 


2 The time is nonAimenslonalited by the time scale. 
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liic core representing the inner region of the spiral initially moves inward and 
upward to conserve impulse. As the core approaches the centroid of vorticity for the 
complete vortex system, - 0.7854, the upperwanl movement ceases and begin 
moving downward while asymptotically approaching x 0 . 

Figure 5-8 shows the amount of circulation contained by the sheet outboard of a 
given panel. The horizontal axis gives the panel index, and since the panels are 
constant in width, this is covalent to the arc length along the sheet Values of 

circulation are negative since the panels’ k were defined to be positive. The y- 
mterccpt shows the total circulation of the system which has been defined to be 1.0 . 
The first curve, the start-up distribution, is elliptical and the strength of the tip core 
appears os the jump from the end of the curve to T - 0. The lessening of the 
gradient indicates the stretching of a section of the sheet. Stretching is extreme in the 
inner-most region of the spiral as indicated by the nearly horizontal inclination of the 
curve while the outer portion is still relatively unaffected by the rollup. Slight 

bumps in the curves therefore shows regions of differential stretching associated with 
the rollup. 

fhe roll-up of a vortex sheet must conserve circulation and momentum. The total 
circulation computed for the numerical model of the vortex sheet as it deform- with 
time shows a negligible fluctuation of about .01% per time step from the expected 
value of 1.00 . The program maintained this value to better than 0.001% accuracy, 
representing only me fluctuations due to numerical round-off errors. Calculated 
centroid of vorticity for a smooth elliptical distribution is 0.7854 . For the numerical 
model of the initial condition used, the centroid of vorticity was found to be 0.7925 . 
During the subsquent steps in time, the centroid value was found to decrease by less 
than 0.05% per time integration cycle, or time step. The steady decrease in centroid 

appears to be connected with the additional panels used to represent the stretching of 
the vortex sheet. 


Betz [10] has shown that the circulation within the spiral, measured from the center 
of the spiral, should approach a constant distribution. Figure 5.9 shows the 
distribution of circulation versus the distance from the position of the roll-up core for 
each time step plotted in Figures 5-7 and 5 8 along with a plot of the approximate 
solution found by Betz. Each curve represents the sheet for increments of 0.100 in 
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nondiincaiiomii time. Since each cycle of time integration coveni 0.005 in 
nondimcntionai time, each curve represent* the result of 20 time steps. .Some of the 
deviation from Metz's result can be explained by the fact that Metz did not use Kaden's 
similarity solution but placed his hypothetical roll up at the end of the sheet. Figure 

>1() Bhows thc Wow u P of th(! region covered by the spiral. In both figures, the area 
over which thc circulation is integrated increases to the right of the plot; thus each 
curve will converges on thc total circulation of the sheet, which in this case is unity, 
toward thc right of thc plot. Although Metz's solution shows the general trends, thc 

numerical solution exhibits a tendency toward higher concentration of circulation inside 
thc roll-up spiral. 


Figure S-At shows thc geometry of thc sheet at elapsed time of 1.00 , The spiral 
planfom generally agrees with those generated by other investigators, the weakness in 
the flat panel method is apparent in uic tightly wound spiral. As thc spiral tightens 
nrar the center, the number of panels become insufficient to adequately resolve the 
spiral. In this case, there are approximately 11 panels describing thc inner-most turn 
of the spiral. Due to the ambiguity in thc radial position caused by spanning flat 
panels across a tight curve, the end of the sheet has crossed over several turns of the 
spiral (Hgurc 5-12). When this happens, the closed section of the vortex sheet stops 
stretching and becomes roughly equivalent to a core of vorticity. 

WhUe sheet crossing represents a violation of the physical nature of vortex sheet, the 
simulation can become unrealistic much earlier. When thc panels forming one turn of 
thc spiral roughly overlap with thc panels forming an adjacent turn, thc geometry of 
the section have been found to stabilize in this configuration. Since the phenomenon is 
produced by the flat geometry of thc panels, this behavior is uncharacteristic of a 
smooth vortex sheet. For this reason, the width of the panels must be reduced in 
order to accurately resolve thc inner regions of the spiral. 


5.2.2. Roll-Up Produced by an Impulsive Movement by a Disk 

The problem studied by Taylor [27] is demonstarted by performing the eUiptically 
loaded wing wake with axisymmetric vortex bands. The motion of the disk is 
modelled by a flat axisymmetric vortex sheet. A rough approximation of the initial 
roll-up of the tip is used (sec Appendix B). The axisymmetric tip core, unlike thc 


two dimensional wire, lias a self-induced velocity which most be determined using the 
conservation of energy. The subsquent roll up of the vortex sheet Into a vortex ring 
proceeds smoothly as shown in the figures. As with the wing wake, the simulation 
begins with 80 panels with the outer 1.25% of the vortex sheet transformed into the 
initial roll-up core with no subsquent core dumping 

Figure 5-13 shows the roll-up of ft vortex sheet produced by an impulsive motion of 

a flat duk. There ate 10 plots, each representing elapsed time of .05, and the initial 
configuration, with the tip core shown as a small circle, shown in this plot. The 
rollup is similar to the wing wake roll-up with one important difference, the spiral 
translation is augmented by the self-induced velocity. The spiral is behaving 
essentially as a ring vortex, Figure 5-14 shows the circulation distribution 
corresponding to the geometries above. The gap between the first curve, showing the 
initial circulation distribution, and the rest of the curves is due to the decrease in 

panel size and a proportional increase in the number of panels imposed by the program 

in the first few time steps. Tins is performed when the distance between the core 
and the end of the sheet is reduced below defined limits, Othewise, all of the 
features of the wing wake case can be observed. 

Centroid of vorticty for an axisymmetric vortex sheet is the radius of a vortex ring 
with the same impulse. The calculated value for the centroid of a continuous vortex 
sheet is x c » 0.8165; the total circulation should be 1. The panel-method introduces a 
error in the centroid to give the computed value of 0.882 for the initial vortex sheet. 
However, the important criteria for the accuracy of the numerical simulation is the 
invariance of these values as the ro. l-up proceeds. The total circulation was found to 
fluctuate by less than 0.0019b while the centroid of vorticity was maintained to about 
0.059b change per time step where each time step is equal to elapsed time of 0.005 . 
In addition, if the roll-up simulation is perfectly accurate, Betz hypothesized that the 
circulation integrated away from the center of roll-up should approach a final 
distribution, A test of the hypothesis is shown by Figure 5-15. The circulation 
distribution in the spiral, which is the region toward the left, appears to converge on 
a relatively steep slope which refeiccts a high concentration of vorticity. This feature 
is more apparent here than in the two-dimensional case because of the faster rate of 
roll-up. Figure 5 16 shows the blow-up of the circulation distributions in the region 
covered by the spiral. 
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Hfiure 517 shows in detail the sheet geometry »t elapsed time of 0.550 . An 
compared to the wing wake roll-up At the same stage of roll-up, the spirals arc looser 
but more turns of the spirals are present. The Invariant tip roll up core Is shown 
with its calculated cross sectional radius which, as can be observed, is very small; it 
represents the tip 1.25% of the initial vortex sheet. Due to the higher rate of inward 
motion in the spiral, the inner-most turn has already reached the point where the 
panels are inadequate to resolve the curvature and a sheet crossing has taken place. 

5.2.3. Roll-Up of a Helicopter Wake 

The evolution of the wake produced by a single rotation of a helicopter rotor was 
simulated by imposing an approximate rotor loading 

rw “ r s* a - (5J) 

on the initially flat axisymmetric vortex sheet This load distribution is shown 
graphically in Figure 5-19. The downwash profile induced by vortex bands, without a 
tip roll-up model, for this distribution is shown in Figure 5-20. The numerical result 
for the velocity at the center of the disk agrees well with the theory, which comes 
out to be 0.392? if the contributions of all of the infloitesmal circular line vortices 
are integrated analytically. 

Figure 5-21 shows the roll-up simulating the evolution of a wake due to one 
rotation of the helicopter rotor at increments of 0.05 in nondimensional time. The tip 
core contains the circulation of the outer 1.259b of the flat vortex sheet and each time 
step represents nondimensional elapsed time of 0.005 . A prominent feature of the 
vortex sheet geometry is the stationary point near the midpoint on the sheet. This 
feature was also found in a similar investigation by John Kantclis [26] using a 
collection of vortex rings to represent the vortex sheet. It is noticed that the 

steepening of the circulation curve, indicating a concentration of vorticity which should 
eventually lead to another roll-up of the vortex sheet, only develops at the center of 
the axisymmetric vortex sheet. Thus, although the vortex sheet contains regions with 
vorticity of opposite signs, the sheet docs not appear to roll-up into two separate 
vortex rings. The circulation plot (Figure 5-22) shows a significant decay of the 
trough in the circulation curve during the first ten or so time steps. This is a 
numerical error due the inaccuracy introduced by the linear interpolation of circulation 
used in the reconstruction of the vortex sheet after each time step. 


” V y 0f °' C ’ tmMm “ "V Invariance of total circulation, 

aqua tn zero, and the impulaa of the vortex tomb Total clrcul.tion and total imnulas 

round to he conmrv* to within 0.001* and .0021* restively during each toe 
gnu on cycle. In addition, the circulation dlatrlhutlon away from the roll-up center 
‘ erected to converge on a „„a. curve. Figure 523 allow. the di.tr, hotion of 
circulation Integra,* .way from the center of the ro,,-„ p f„ r each time atop to 

•gum, 5-21 and 5-22 overlapp* to .how the comparative deviation* The blow-up of 
the region covered by the spiral is ihown in Figure 5-24. 

The geomc.ry of the vortex .hee, at elapsed time of 050 i. toown to Fig,. re 5-25 

r'r the ”' i,,im " m 10 ** dre “ kti '”> dtotributlon curve (Figure 

Cire “ “ thow «“ “““ Of «* negative vorticity to to. roll-up. Tip 

* 1. much facer than the two previous cm« to.ee toer, fe mom circulation 

concentrate! to the tip region. The blow-up Figum 5-26. tone, ertenaive tom, ousting 

tor the inner-most turns of the spiral. Tin cor* 

.... _ v p lip coie » whlch “PJPcars as a dot in the 

nnddlc of the spiral, is shown with its cross-sectional radius in scale. 

3.2.4. Kclvin-Hclmholtz Instability 

WUdl * ,he of underlying Kalvto-Halmholtx 

bdtty becomes evident under amend condition* A wake vonex tom, which 1, 

Z” finely discretized touched toe geometry ahown to Figure 527 after ten time toep* 

Dunng each thue atop. to. juneto werc .How* to move to, maxtoium of 0.1 time. 

c panel width. The circle on the tight reprcsents to. location of Radon's corn model. 

“ me * M ta "" Ver,iC,a dlrK, ‘° n> *• <leVdt ^“* of Kelvto- 
mholta-type toStobUlty becomes visible (Figure 5-28). Further magi„ifi«i„n .how, 

the ™v« to have a period of „ panto, M d toe ecelto * be witoto to. » digit 
accuracy of the computer. ^ 

The damping „f Relvta-Holmholtz Instability to toe numeric^ eimulatlon i, m toe 
orm o ow-pasa filtering the waving in the vortex tomt, Thi. filtering occnia aa 
panels apan KroB crest, and , tough, of .hot, wavelength wave* Due to the 
atretobmg of toe torn,, the panto, which describe geometry of toe vortex tom, 
following a time integration cycle win he .lightly di.pl.eed from toe previou. panel* 
amount of torn displacement determines the amount „ f m « rln g which take, place 
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This is becauar wavincss in the displacement of ft ftheet formed by panels can not be 
resolved unless the the next net of panels overlap the displaced panels, For this 
reason, the instability appears when the amount of displacement, limited by the 
amount of time elapsed in one cycle, is kepi, small. This is what happened in Figure 
5-28. 

The filtering effect can also be suppressed by continuously increasing the spatial 
resolution of the panel-method. An example of this numerical behavior is shown in 
Figure 5-29. The initial conditions on the sheet and the time integration cycles arc 
identical to the disk wake roll-up case presented earlier. To illustrate the filtering 
effect, the size of the panels were constantly reduced to about one half that of the 
panel used for Figure 5-13. This is accomplished by specifying GFAC ■» 2 instead of 
GFAC - 1; see Appendix C. This effectively suppressed the filtering and the small 
instabilities are able to grow more freely. As the roll-up proceeds, the secondary roll- 
up spirals are amalgamated into the tip roll-up spiral and the secondary spirals are 
also stretched as they spiral into the tip roll-up spiral. This behavior can be observed 
directly above the main spiral where stretching has deformed the secondary spiral. 
Plot of the circulation distribution, Figure 5-30, shows jumps in circulation which 
corresponding to each secondary roll-up. The circulation distribution in this plot was 
initially elliptical like those in Figure 5-14. Except for the existence of the secondary 
roll-up, the simulation proceeded no differently than the previous case. The computer 
program maintained the invariance of circulation and momentum with the same degree 
of accuracy as before. Distributions of circulation taken from the center of the main 
roll-up spiral is shown in Figure 5-31. Figure 5-32 shows the blow-up of the region 
covered by the main roll-up spiral They are no different from the previous results 
except for the effects of secondary roll-up in the circulation distribution. 
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Start: 

Program 'TOI.ICAU ." 


Initialize: 


*1J 

Panel Position 


Panel Orientation 


anti Width 

Tij 

Panel Strength 

( i : 

Panel Index, j : Time Index) 


"JNITIA 


Calculate (x i0 , Dx 10 , 7 10 ) 


* Subroutine names 
arc in quotes. 

• ’RUNGI'l” uses the 
subroutine "VJfljQCX’ 
wliich calls the 

velocity routines: ’'VJiI.21) 1 ’, 

w v}iUN ,, ,"viiLOu , r, M vj-:uo 


Data Output: 

Hie #1 Cj» imax^ 

Hie #2 (xjj, Dx^j) 

File #3 (i, *v< j) 


Time Integration: 

j+i m fxnCx^j, \y% j j , 7t j) 


j - j + 1 


Segmentation: 

ki + y 2 j+l, 7ij + l) - fxnCxij+2, 7ij) 


^ j " inax 
then STOP. 


Sheet Rcdiacretkation: 

+ “ fxn(xj + y 2 j +1 ) 


Figure 5-1: Flow Chart of the Vortex Dynamics Program 
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6. Combine Velocities, 

eg- v ic * (vij + Vij + y 2 ), for 
Translation to the Final Position.: 
*1J+1 " *ij + v ic At - 



Figure 5-2: The First Order Kunge Kutta Scheme 


in fit a spline Ijctween 
points i and ill. 

1. Obtain adjacent jxtints 
(Extrapolate at the ends). 


2. Slopes at i and il l 
arc obtained front the 
tangents to the arcs 
G 'Uitl) and (Mt l»it2). 


3. A parabola is fitted 
between i and i+1 
to fit the dopes. 




S • 
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Figure 5-3: The Reconstruction of the Vortex Sheet 


Subroutine "RUNGEl” gives the 
translated panel centers, x 1 j + 1 . 
Subroutine "REDISK" splines the 
panel centers and determines the 
spatially equidistant points, 
along the spline. 

Subroutine "PANELKR” spans points 
x f"VitJ+i *i+y 2 j + i with panel i. 




Figure 5 4: lire Rcdiscrptization of the Vortex Sheet 
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I’X|U»1 width panel discretization 

t»r a flat wake vortex shrct. ” -» — * — ■► 


*l'he remaining sheet is 
stretched (rediscretized) 

so the edge is under the core. — "* — 


KtJ*« Solrt , 

f- i i t 


$ 


Figure 5-5: Initial Configuration of the Vortex Sheet and Roll-Up 
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Figure 5-29: Repeat of Fig. 5~13 with Finer Paneling: T « C 



Figure 5-30: The Corresponding Circulation Distribution 
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Chapter 6 
Conclusions 


The behavior of vortex sheets with axisymmetric geometry presents the next frontier 
in the study of vortical flows. Asymmetric flow Ues between the relatively well 
understood two-dimensional flows and the complex three-dimensional tlows. For this 
reason, understanding of axlsymmetric vortex sheets is useful before completely three- 
dimensional vortex sheets can lx modelled. Due to the complexity imposed by the 
increase in number of variable as the number cf dimensions arc increased, three- 
dimensional models using discrete vortices or dipoles are more desirable than panel 
methods. However, the use of simpler models would be meaningless unless the 
accuracy of the models can be determined. For this reason, the axisymmetric panel 
method was developed as a higher-ode r model to test the accuracy of simpler 
axisymmetric vortex models. 

To accurately formulate the behavior of an axisymmetric vortex panel, a flat strip of 
a uniform two-dimensional vortex sheet, called the vortex panel, was first studied. 
Tin calculations of the velocity field show that an assemblage of panels will mutually 
induce motions which are more representative of a vortex sheet than an assemblage of 
point vortices. Even when reproducing a vortex sheet with non-uniform vorticity and 
curved geometry, the panel method performs no worse than the discrete vortex method. 
The simple geometry of the vortex panel allows a straightforward derivation of a 
matched asymptotic solution for the flow induced by an axisymmetric equivalent ol 
the vortex panel, the vortex band. Vortex bands can be assembled to accurately 
represent an axisymmetric vortex sheet, limited by the geometric inaccuracy in 
resenting a smooth curve using straight line segments and, at least for this study, 
the non-smoothness of the vorticity distribution. Both of these weakness are problems 
of spatial resolution which can be alleviated by increasing the number of elements. 
Th USl panel-method can be used to accurately simulate the evolution of an 
axisymmetric vortex sheet. 
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Die voitex sheet dynamic* calculation is an initial condition problem. One Important 
rc«m« lma been the accuracy nf the panel- method lu rep.iKb.dng the llnw induced by 
« flat vortex sheet which is the initial state ol' the vortex sheet in many problem of 
interest. In the present atiuly, the singularity at the edge of the vortex ahcct won 
modelled using a simple discrete vortex. While the introduct* of a tip vortex core 
waa necessary for the accurate reproduction of the glottal flow field, it is also useful 
ior the numerical simulation of the self Induced deformation of the vortex sheet. 

The self induced velocity of curved vortices have inhibited the atudy of the dynamic 
of axlsymmetric vortex sheets. My using the method of matched asymptotic expansion, 
the self Induced velocity of an axisymmetric vortex panel consistent with the kinetic 
energy of a section of an axisymmetric vortex sheet was found. The axisymmetric 
panel-method was u.scd to simulate the roll up caused by the impulsive motion of a 
circular disk and the roll-up of a vortex sheet with a loading which approximates that 
of a helicopter rotor. In both cases, smooth roll-up was observed and both total 
circulation and impulse were conserved. My relaxing the lovrfw, filtering inherent in 
the numerical method, secondary roll-up due to Kelvin- Helmholtz Instability can be 

observed. This feature for an asymmetric vortex sheet was simulated numerically 
using the panel-method described above. 


78 


References 


Rcwenhcad, !.«, w l'hc formation of Vortices from a Surface of Discontinuity'*, 
ProcJiSocJxtndon Ser.A, Vol. 134(1931), pp. 17<F92. 

Eirkhoff, G n Fisher, J., "Do Vortex Sheets Roll Up?”, Clrc»Mat. Paler no, Vol. 
8(1959), pp. 77-90. 

Saffman, P.G., Baker, G.R., "Vortex Interactions”, Ann. Rev. Fluid MecK Vol. 
11(1979), pp. 95-122. 

Moore, D.W., "On the Point Vortex Method”, SIAM J£clStaX.Comp n Vol. 
2(198”), pp. 65-84. 

Baker, G.U., The Cloud-in-Cell Technique Applied to the Roll-Up of Vortex 
Sheets", J.Comp.Physlcs, Vol. 31(1979), pp. 76-95. 

Murman, E.M., Strcmel, “A Vortex-Wake Capturing Method, for Potential Flow”, 
AJAA Paper, No. 82-0947, iime (1982). 

Moore, D.W., "Ihe Spontaneous Appearance of a Singularity in the Shape of an 
Evolving Vortex Sheet”, ProcJMocAondon SerA, Vol. 365(1979), pp. 105-119. 

Meiron, D.I., Baker, G, Orzsag, S.A., The Analytical Structure of Vortex-Sheet 
Dynamics. Part 1. Kelvin-Ilclmholtz Instability”, J Fluid MecK Vol. 114(1979), 
pp. 283-298. 

Karman, Th. von, Burgers, J.M., Aerodynamic Theory, Springer, (1935), pp. 
320-325, See also Westwater 12 

Betz, A, "Uehav’or of Vortex Systems”, NACA TM~713il933). 

Pullin, D.I., "On a Generalization of Kaden’s Problem”, /. Fluid MecK Vol. 
104(1981), pp. 45 53. 

Westwater, 1.I-, "Rolling Up of the Surface of Discontinuity Behind an Aerofoil 
of Finite Span”, AeroMes.CommJlAM^ No. 1692(1935). 

Fink, P.1,, Soli, W.K., ‘^Calculation of Vortex Sheets in Unsteady Flow and 
Applications in Ship Hydrodynamics”, Proc. 10th Symp. Naval Hydrodynamics, 
Cambridge, Mass., (1974), pp. 463-488. 

Fink, P.T, Soh, W.K, "A New Approach to Roll-Up Calculations of Vortex 
Sheets”, Proc.RSixJ/ondan Ser.A, Vol. 362(1978), pp. 195 209. 

Baker, G., A Test of the Method of Fink & Soh for Following Vortex Sheet 
Motion”, /. Fluid MecK VoJ. 100(m;s). pp. 209-220. 


79 


16. Ilndjmukers, H.W.M., Vaatstta, "A Higher Order l'and Method Applied to Vortex 
Sheet Hoi 1 VJp", Al AA Fajw, No. 82 0096(1982). 

17. Higdon, J.L., Po/rikidls, C„ The Self-Induced Motion of Vortex Sheets” /. Fluid 
Med i, Vol. 150(1985), pp. 203 231. 

18. Murky, M., “Calculation of Vortex Sheet Roll-Up in a Rectangular Wind 

Tunnel", / .Aircraft , Vol, 12, No. 9(1975). 

19. Batchelor, U.K., Introduction to Fluid Dynamics, Cambridge University, (1957). 

20. Widnalt, S.H., Tsni, C.Y., The Instability of the Thin Vortex Ring of Constant 
Vortieity”, PhtLTran.PSocJxmd^ Vol. 287(1977), pp. 273-305. 

21. Lamb, II, Hydrodynamics, Dover, (1932). 

22. Bliss, D, The Dynninius of Curved Rotational Vortex Lines”, Master’s thesis, 

M.I.T, Aero/Astro, (1970). 

23. Bender, C.M., Orsssag, S.A., Advanced Mathematical Methods for Scientists and 
Engineers, McGraw-Hill, (1978). 

24. Jabnke, E, Emde, I ; , 7 'able of Functions, Dover, (1945). 

25. Tung, C, Ting, L, “Motion and Decay of a Vortex Ring”, The Physics of 

Fluids, Vol. 10, No. 5(1967). 

26. Kantelis, J.. Private Communications 

27. Tayor, G.I, "Formation of a Vortex Ring Giving an Impulse to a Circular Disk 
and then Dissolving it Away”, J. Appl. Fhy, Vol. 24(1953), pp. 104. 

28. Roache, PJ, Computational Fluid Dynamics, Hermosa Publishing, (1972). 


80 


Appendix 


The following pages are divided into three sections, 

1) The listings of the FORTRAN routines for competing the velocity induced 
by various vortex elements. 

2) The description of the numerical model of the initial state of the vortex 
sheet prior to the numerical roll-up simulation and the FORTRAN listing. 

3) The listing of the FORTRAN program for the numerical simulation of vortex 
slret dynamic, 

4) The inputs for the vortex sheet dynamics simulation program used to 
generate the figures within the text. 


All of the following program listings have been written in FORTRAN 77 
Version 3.0 for the Digital Equipment Corporation VAX-750 computer. 
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Appendix A 

The Vortex-Velocity Routines 

The velocity indeed by a vortex element is calculated uving a Cartesian coordinate 
system, (X,Y), centerH on the vorte v inducing the flow. The panel’s inclination 
respect to the coordinate system is expressed by (DX,DY), where DY/DX equals the 
tangent of the angle between the panel panel and the x-axis. (W) is the half-width 
of the panel and (G) is the vorticity along the panel. The velocity is obtained as 
orthogonal components, (u,v). 

The names of the two-dimensional velocity routines: 

VEL2D : The Vortex Panel Routine 
VEL2DF: The Point Vortex Routine 

The axisymmetric vortex routines differ from the two-dimensional vortex routines due 
to the inclusion of the radius of the vortex element, R. The composite solution, V cornj)1 
for the velocity induced by a vortex band is 

V - V + V - V 
comp in T out Mo' 

where the routine 

VELIN computes the inner-solution, V ln 
VELOUT computes the outer-solution, V QUt 
VELIO computes the intermediate-solution, V lQ . 

(Note: The V out alone gives the velocity induced by a simple vortex ring.) 


SUBROUTINE VEL2P (UIN, VIN) 


* THIS SUBROUTINE CALCUI.ATES THE VEl.OCITr FIELD INDUCED 

* BY A FLAT VORTEX STRIP OF 

* WIDTH : 2W 

* STRENGTH : C 

* INCLINATION : ( DX, DY ) 

* AT POSITION <X, Y) RELATIVE TO THE STRIP MIDPOINT 

* (DUMMY VARIABLE: RADIUS : R) 

C RELATIVE COORDINATES AND QUADRATIC-POSITION VARIABLES 

C 

CC*MQN A£LO/ R. X, Y, DX, DY. W, G 
PI - 2«ASIN( i . ) 

COST - DX / SQRT(DX..2+DY**2) 

SINT - DY / SQRT(DX**2+DY**2) 

05 - (X*»2 + Y**2)/W»*2 
Q4 ■ <Y*COST - X«SINT)/W 
03 • (X*CQST + Y*SINT)/W 
02 » 05 - 2.Q3 + 1. 

01 » Q5 + 2*03 + 1. 

CHECK FOR SINGULAR CONDITION 
IF (Q4.EQ.0) THEN 
ATN - 0.0 
CO TO 740 
END IF 


ATN * ATAN((1 .=G3)/Q4)+ATAN((1 .-K)3)/Q4) 
740 CONTINUE 

IF ((01 . EO.0) .OR. (Q2.LQ.0)) THEN 
UIN - 0 
VIN - 0 
CO TO 750 
END IF 

C PERTURBED VORTEX STRIP FORMULA 

U1 —SINT«L0G(02/Q1)42.*COST*ATN 
VI •-COST*LOC(Q2/Q1)-2.*SINT*ATN 
UIN *> ~C / 4./ F>I • U1 
VIN » C / 4./ PI • VI 
750 RETURN 



« o 
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SUBROUTINE VE' 2DF (U1N. VIN) 


THIS SUBROUTINE CALCULATES THE VELOCITY FIELD INOUCPO 
BY A FLAT VORTEX STRIP OF 

WIDTH : 2 W 

STRENGiH : 0 

INCLINATION ; ( OX, DY ) 

AT POSITION (X. Y) RELATIVE TO THE STRIP MIDPOINT 
(DUMMY TRIABLE: RADIUS : R) 


RELATIVE COORDINATES AND QUADRATIC-POSITION VARIABLES 

COMMON /VELO/ R. X, Y, DX, DY, W, G 

PI » 3.141592654 

05 - SQRT(X**2. + Y«*2.) 

IF (Q3.EQ.0) GO TO 730 
COST * X / 05 
SINT - Y / 05 


U1 - SINT/Q5 
VI - C0ST/Q5 
UIN —0/2. /PI • U1 
VIN - C/2. /PI • VI 
750 RETURN 
END 
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SUBROUTINE VELIN (USN, VIN) 


THIS SUBROUTINE CALCULATES THE INNER SOLUTION 
OF THE VELOCITY FIELD INDUCED BY AX I SYMMETRICALLY 
PERTURBING A FLAT VORTEX STRIP OF 
WIDTH : 2W 

STRENGTH : 0 

INCLINATION : ( DX, DY ) 

RADIUS : R 

AT POSITION (X, Y) RELATIVE TO THE STRIP MIDPOINT 

THE DERIVATIVES OF THE STREAM FUNCTION HAS BEEN 
DIVIDED INTO TWO PARTS: 

1 - RECTILINEAR SOLUTION 

2 - FIRST ORDER PARTICULAR SOLUTION 


C RELATIVE COORDINATES AND QUADRATIC-POSITION VARIABLES 

C 

COMMON AELO/ R, X, Y. DX, OY, W, G 
PI - 3.141592654 

COST - OX / SQRT(DX*«2.+DY**2.) 

SINT ■ DY / SQRT(DX*«2.+DY*»2.) 

Q5 m (X*»2. + Y»*2 . )/W»*2. 

Q4 - (Y*COST - X«SINT)/W 
Q3 - (X*COST + Y»SINT)/W 
02 - Q5 - 2.. 03 + 1. 

Q1 - 05 + 2. *Q3 4 1. 


CHECK FOR SINGULAR CONDITION 
IF (Q4.IQ.0) THEN 
ATN « 0.0 
GO TO 740 
END IF 


ATN » ATAN((1.-03)/Q4)+ATAN((1.403)/Q4) 
740 CONTINUE 


C PERTURBED VORTEX STRIP FORMULA 

U1 — SINTeL0G(Q2/Q1)+2.*CQST»ATN 
VI — C0ST»L0G(Q2/Q1)-2.*SINT*ATN 
U2 - (X/2./R) o U1 
V2 - (X/2./R) • VI 

• + .5*W/R •( L0G(Q2*Q1 ) 

• -03 • LOG (02/01) 

• +2.*Q4»ATN 

• -4.) 

UIN — C / 4. / PI . (U1 + U2) 

VIN * G / 4. /PI • (VI + V2) 

750 RETURN 
END 
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SUBROUTINE VELOUT (IK), VO) I ELLIPTIC INTEGRALS i 


THIS SUBROUTINE CALCULATES THE VELOCITY FIELD 
INDUCED BY A VORTEX RING WHICH *5 USED AS THE 
OUTER VELOCITY SOLUTION FOR A VORTEX BAND OF 
WIDTH : 2W 

STREJGTH : 

INCLINATION : (DX. DY) 

RADIUS : R 

AT POSITION (X. Y) RELATIVE TO THE BAND MIDPOINT 


COMMON /VEL 0/ R, X, Y. DX, DY, W. C 

REAL Ml. LAM, K0, K1 . K2, KL0, KL1 , KL2, K, KP 

uo * t.e 

vo « 0.0 

PI » 3. 141592654 


R1 *• SQRT(X**2.+Y**2. ) 

R2 «* SQftT((X+2.*R)**2.+Y**2.) 

LAM - (R2-R1 )/(R1+R2) 

Ml • 1. ~ LAM* *2. 

IF (R1.EQ.0) GO TO 730 

IF (LAM.EQ.0) GO TO 730 

DXR1 - X/R1 

DYR1 - Y/R1 

DXR2 *■ (X+2. *R)/R2 

DYR2 - Y/R2 

DXL - (DXR2-DXR1 )/(R2+R1 ) - LAM/(R2+R1 )*(DXR1+0XR2) 
DYL « (DYR2~DYR1)/(R2+R1) - LAM/(R2+Rl)*(DYR1+DYR2) 
K0- 1.3862944 
K1 - .1119723 • Ml 
El - .4630151 • Ml 
K2 - .0725296 • Ml*. 2. 

£2 - .1077812 * M1**2. 

KL0 - .5 

KL1 • .1213478 • Ml 

ELI - .2452727 • Ml 

KL2 •> .0268729 * Ml*. 2. 

EL2 - .0412496 * Ml* *2 

K - K0 + K1 + K2 - (KL0 4 KL1 + KL2)*L0G(Ml) 

E - 1. + El + E2 - (ELI + EL2)*L0G(M1) 

KP • £ / LAM / Ml — K / LAM 

EP «(E - K) / LAM 

DXR - 0XR1 + DXR2 

OYR - DYR1 + 0YR2 

RR » R1 + R2 


UO - C*W/PI/(R+X) • (DYR* (K-E)+RR* (KP-EP) *DYL) 
VO -~C*W/PI/(R+X) • (DXR* (K-EJ+RR* (KP-IP) *DXL) 
RETURN 
END 


730 


fl".OT«rjwai«Trin*" 
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SUBROUTINE VELIO (UIO, VIO) 


THIS SUBROUTINE CALCULATES THE OUTER LIMIT 
APPROXIMATION Or THE INNER VELOCITY SOLUTION 
GIVEN BY ROUTINE VEEIN IN OR ER TO MATCH THE 
OUTER VELOCITY SOLUTION GIVEN BY ROUTINE VCLQUT 


COMMON /VELO/ R, X, Y. OX, DY, W, G 
PI * 2* A5IN(1.) 

COST «■* PX / 5QRT(PX**2+DY**2) 

SINT « DY / SQRT(&X**2+DY**2) 

Q3 *» (R**2 > Y»*2;/W*»2 
Q4 « (Y*COST + R*SINT)/W 
03 - (-R«COST + Y*S1NT)/W 
02 « 05 - 2*Q3 + 1 . 

Q1 - 05 + 2*03 + 1 , 

UIO - 0.0 
VIO - 0.0 

R1 - SQRT(X«*2.+Y**2.) 

IF (R1.EQ.0.0) GO TO 720 


UIO —G*W/PI • ( 1. + X/R/2 . ) • Y/R1/R1 

VIO * G*W/P I • (( 1. + X/R/2.) • X/R1/R1 + .5/R*L0G(R1/8/R)) 

720 VIO - G*W/PI • UOG( 8*R/W ) /2 /R 
• + VIO 


RETURN 
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Appendix B 

The Initial Conditions for Roll-Up Simulations 

The velocity field induced by ft distribution of panels alone precludes ft nmooth toll 
up unless the singularity at the edge of the vortex ahect is correctly modelled. The 
vorticity singularity at the edge of a vortex sheet is expected to produce an 
infinitcsmal roll-up in an mfinitesmal amount of time. Such a roll-up relieves the 
singular velocity discontinuity which would otherwise exist at the edge of the vortex 
sheet. 

The roll-up simulations uses an initu configuration which represents a vortex sheet 
with a model for an arbitrary amount of roll-up at the edge of the sheet. The 

initially rolled-up section of the vortex sheet is replaced by a circular core of vorticity. 

For the roll-up of the wake vortex sheet shed by an clliptically loaded wing, Kaden’s 
similarity solution for roll-up of parabolic loading is used to approximate the position 
of the center of the roll-up spiral. The position of the center of the roll-up spiral is 
given in Reference [9J. Figure H-l shows Kaden’s solution for the roll-up of the outer 
X 0 segment of the vortex sheet. 

a 0 ■ 0.S7 X 0 
b 0 - 0.88 X 0 

A significant feature of Kaden's solution is the stretching of the remainder of the 

sheet up to the point directly under the center of the spiral, as noted by the 

displacement of point A on the sheet to point A’. This displacement conserves the 
overall centroid of vorticity since the center of the spiral is less tuan the centroid for 
the outer X Q section of the flirt vortex sheet. The numerical model replaces the 
rolled-up X 0 section with a core at the center and linearly stretches the remaining flat 
vortex sheet so its edge would come directly under the core. The resulting model 
reproduces the centroid of vorticity for the elliptic loading, x c « tr/4, to within 1% 
accuracy. 


> .* 
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While Kaden's method In available to approximate the location o' a tip core for two 
dimensional roll up, its application to an axlsyminetrk roll up will not conserve 
impulse, lherefore, a vortex ring with the sume impulse and energy art the outer X f 
aeetlou ol the axisymmetric vortex sheet is used oh a core* However, the vertical 
location ol the rollup core can not be obtained by the conservation equations. Since 
no theory exists for obtaining the vertical location, the ratio of distances for Kadcn’s 
solution was used to approximate the location. In addition, since the core conserves the 
centroid of vorticity, the remuini, g sheet was not stretched for this cose. Figure H 2 
shows the numerical model of the initial roll-up for ax I symmetric vortex sheet which 
represents a translating disk, 


~ X t 1 

C. - J, l R 


r « SK exp l 


2X X 2 .. 
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4(2KX-Xl) V * 4 

where, R : the radius of the original axisymmetric vortex sheet, 
r Q : the radius of the roll up core. 


The initialization of the simulation of vortex sheet dynamics is performed by the 
routine INITiA. 


Hguro B-l: Kndcn’s Similarity Solution and its Numerical Model 



Figure R2. Mtial Roll-Up Model for Axisymmctric Vortex Sheet 
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SUBROUTINE INITIA (NCOfil ) 

INITIAL VORTEX SULLY GEOMETRY FOH PROGRAM TULLCAGl. 

modified for complete Kuden'e ■oImUoii ond OKleymmetrlc tip florae. 

Common block "PANEL" give* th# poeltlon* (X,Y), the loci loot Ion (DX.UY), 
the width (W) and th# gradient of clrcw'atlon (DG), 
Common block "TIME" give# number of panel* at eauh time etep and the 
time nlup data. 

Common block "CORE" giver, the Information regarding the tip roll -ap core. 
Common block "GAMMA” give* the poeltlon of each panel on th* elliptical 
circulation plot (GS.G). 

Common block "ZOT” epeclfli# th* typ* of roll-up: 0 « 21) elliptic wing 

i 1 m uniform dl*k wok* 

I 2 « model rotor wake 

COMMON /PANEL/ X(0: 505), Y(0:505), DX(0:505), DY(0:505), 

# Vf(6 ; 505) i 00(0:505) 

COMMON /TIME/ NPANEL, NSTEP, NTIME, TIME, TSTEP 
COMMON /CONE/ XC, YC, GO. RC, UC, VC 
COMMON /GAMMA/ 0(0:505), GS(0:505) 

COMMON /ZOT/ ZOT 

PI - 2.*ASIN(1 .) 

X(0) - 0 
Y(0) * 0 
G(0) - -1. 

IF (ZOT. OT. 1) G(0) ■» 0.0 
GS(0) » 0 

I Equal width eegmentat I on routine. 

I R0 and Rt ore the poaltion* In span of the Inner and the outer edge* of 
I the panel respectively . 

DO 10 I ■ 1, NPANEL 

X(I) - (FLOAT(I) - .5) / NPANEL 

Y(I ) » 0.0 

R1 - FLOAT(I) / NPANEL 
R0 * (FLOAT(I) ~ 1.) / NPANEL 
IF (ZOT . NE . 0 . AND . ZOT . NE . 1 ) THEN 
SF (R0.EQ.0) THEN 

DG(I) - -(Rl»*ZOT) • SQRT(1 .=R1**2)*NAPNEL 
ELSE 

DG(I) - -"(Rl »*ZOT*SQRT( 1 ,~R1**2) - R0*.ZOT*SORT(1 ."R0**2))»NAPNEL 
END IF 

G(I) - -P1**Z0T*SQRT(i .-R1..2) 

ELSE 

DG( I ) - ~(SQRT(1.~R1*»2) - SORT(1 .-R0»«2))»NPANEL 
0(1) « -SQRT( 1,- R1 *»2) 

END IF 

GS(I) " X(l) 
ff( 1 ) m .5/NPANEL 
DX(I) « 1. /NPANEL 
10 0Y( I ) m 0.0 


I Tip Core Routine* 

I 2D *lllptlcally loaded wing wake. 
IF (ZOT.EU.0) THEN 


-v 



xc - i. - .&■/ • ( float(noorf) / npanli ) 
yc - .aa * (, float (ncoro / npanli. ) 
cc •» -c(npani:l~ncorl) 

PKC - KC - X(Nf'ANI |,~NCOHi:) - W(NPANFI.' NCWT) 

W J 5# 1 " 1 . npani:l ncorl 

X(l) - DXC * FLOAT (M)/(NPANi:E NCORE) 4 DXC/VtW'ANM.-NCOliC) 4 X(n 
W(I) " DXC /(NPANLI' tit ONE) 4 W(|f) 

DG(|) M 1,0(1) - C(I 1)) /i Ml) 

SO CONTINUE 

I 

I Impulsively moved fiot dleh wcke 


ELSE IF (ZQT.UE.1) THEN 

WC « FLQAT<NCQHE) / NPANLL 

XC - ?QRT< 1. - 2..WC/3. 4 (WC.*2)/3.) 

YC « .88 / .57 * (1. - XI.) I This l« * guess; see Kaden. 

CC « -G(NPAhEL-NCQRC) 

( RC - a. • XC / EXP( PI. *2 /4 .XC /SORT<?. .RC - WC..2) + 7./4.) 

I Model of a helicopter rotor wake 


ELSE IF (20T.EQ.2) THEN 


I 

t 

I 

I 


tlC m FLOAT (NCORE) / NPANEL 
XC » SQRT( 1. - 2..IVC/3. 4 «C.*2/3.) / (1 
YC * ,8B / .37 • (1. - XQ) 

GC « -G ( NPAN E L“<KCORE) 

RC « 8. »XC /iiXP( PI..2 /4 *R /(1.~KC)*.2 /SQPT(2 ,.WC -#C*»2) 47. /4.) 


■WC) 

Thle 


1* only a gu*g«, 


END IF 

I Initial Condition Record 


30 


WRITER,.) 'INITIAL SHEET: I. X(I), Y(I), W(I), DG( I) * 

DO 30 I - 0. NPANEL 

GS(I) - SQRT(1. - C(I)*.2> 

*RITE(8.») I, X(l), Y(l) , Wi I > . DC(I) 


WRITE(Q,*) ‘CORE: 
WHITE(8,») • * 


NCORE, XC. YC. GC. RC’, NCORE, XC. YC. GC. RC 


NPANEL - NPANEL - NCORE 


RETURN 

FND 
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Appendix C 

The Vortex Sheet Dynamics Program 

The program TOIXCAU. combines the subroutines listed below to simulate the self- 
induced motions of a vortex sheet. Input file for the program is stored as device #7. 
The input file is in the form: (RN, NPANEL, NSTBP, TFAC, GFAC, NCORE, ZOT); RN 
is the run’s identification number, NPANEL is the number of panels initially used to 
describe the vortex sheet, and NS'Ilil* is the number of time integration cycles for the 
run. TFAC defines the length of time used in the time integration, TFAC * 1 is 
equivalent to the nondimensional time span of .005. GFAC defines the size of the 
panel relative to the distance between the sheet and the center of the roll-up spiral. 
It is used to insure the accuracy of the tip region as the distance between the sheet 
and the roll-up core is diminished during the roU-up. NCORE specifies the number of 
panels describing the initial vortex sheet to be amalgamated into a tip roll-up core. 
The type of rc.’J-up, as determined by the selected initial condition, is controlled by the 
parameter ZOT. 

ZQT - 0: Two-dimensional roU-up of wake vortex sheet shed by an 
elliptically loaded wing. 

ZOT - 1: The roll-up of a vortex sheet simulating the formation of 
a vortex ring by the impulsive motion of a circular disk. 

ZOT - 2: Axisymmctric roll-up of a circular vertex sheet model 
of the v/ake vortex sheet shed by a hclicoptov rotor. 

The program version presented here outputs two files describing the configuration of 

the vortex sheet after every ten time integration cylces. The file containing the 

geometrical data of the sheet is in devi;*. #9 and the file containing the distribution of 

vorticity is in device #10. Device #8 holds a utility file describing the execution of 

the program. 


i 
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PROGRAM TOEECAEl 


EXECUTIVE ROUTINE FOR ROEEUP PROGRAM 

COMMON /PANEL/ X(0:305), Y<0:505), DX(0:505), DY(0:905), 

J W(0:S05), DG(0:5B5) 

COMMON /CONV/ U(0: 90S). V(0: 305) 

COMMON /TIME/ NPANEE, NSTEP. NT1ME, TIME. TSTEP 

COMMON /GAMMA/ 0(0:505), CS(0:5Co) CS ■ « not 

COMMON /CORE/ XC. YC. GC, RC, UC. VC I *" tM * v«rtlon. 

COMMON /TFAC/ TFAC, UFAC 

COMMON /ZOT/ ZOT 


Input data read: 


a. Reference Number (IRON) 

b. Number of pan* It initially In th# thttt (NPANEE) 
t:. Maximum numbtr of time tttpt In tht run (NSTEP) 
d. ( ,005/TFAC) !• tht els* of toch tlm« tltp 

t. (J.h.et edge to tip cortl/GFAC) It tht pontl width 
f. Number of cortt Initially dumped Into the tip core 
a. Type of loading (ZOT) 


READ (7 , •) IRUN, NPANEE. NSTEP, TFAC, GFAC, NCORE. ZOT 
NINIT • NPANEE 


Initial Izt output f I It 

WRITE(8, •) ’ROLL UP PROGRAM FOR VORTEX SHEET. 

WRITE(8.*) ’RUN NIM3ER: IRUN 

WRITE(8,*) 'SPLINING BY ROUTINE REDISK’ 

WRITE(B.e) ’INTEGRATING TIME FACTOR : ’. TFAC 
WRITE(B,.) ’GEOMETRIC WIDTH FACTOR : *. GFAC 
WRITE(B.«) ’INPUTS: (NPANEE. NSTEP) *. NPANEL, NSTEP 
WRITE(8,*) ’ZOT FACTOR: LOADING TERM. AXISYMrfETRIC IF >1 
WRITE(B.*) * ’ 


ICHIRO SUGIOKA 37-481 ’ 


ZOT 


Compute Initial Condition 
CALL INITIA (NCORE) 


I — - - ■ — — — 

I Time Integration Loop 

DO 10 NTIME « 0. NSTEP 

CALL WRITER l Write output data file 

CALL RUNGE1 I Rlnge-Kutto time Intgratlon tcheme 
I CALL EULER I Euler time Integration tcheme option 

DELTAL - 2. • W(2) I Panel width parameter (OELTAL) 

I m to a default panel width 

CALL REDISK (X.Y.G, NPANEL, NPANEL, GMAX.SUML, DELTAL) 

I Spline* ond dlvldee the theet 

CALL PANELER(X,Y.DX.DY,W.DG,C, NPANEL) 

I Forme the panel* 


10 CONTINUE 

_ , End of Loop ........ — 

CALL WRITER I Record tht end retulte 

END 
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SUBROUTINE WHITER 

I Thin program record* th* r**ult* of tollcall ot> fllu 
l FOR80B.DAT (Control accounting Information) 

I FOR008.DAT (Geometrical grapMc* data) 

I FORO10.DAT (Circulation graphic* data) 

COMMON /PANEI / X(0:5d5) ,Y(0:505) ,DX(0:505) ,OY(0:505) . 

| W(0:50S), 00(0:505) 

COMMON /TIME/ NPANEL, NSTEP. NTIME, TIME. TSTEP 
COMMON /CAMMA/ K®:»05). CS(0:505) 

COMMON /CONV/ U(0:505). V(0:505) 

COMMON /CORE/ XC . YC. GC, RC. UC, VC 

WRITE(B,*) ‘RESULTS FOR NTIME. TIME, NPANEL: * , NTIME, TIME, NPANEl 

IF (FLOAT(NTIME/10).NE.FLOAT(NTIME)/10.) CO TO 800 

WRITE(9,») NTIME. NPANEL 
WRITE(10,«) NTIME. NPANEL 

i Writ* g«om*trlcal data for lat«r graphic*. 

DO 10 I - 0. NPANEL 

WRITE(9.*) X(I), Y(I). DX(I), DY(I). W(I) 

10 CONTINUE 

WIiITE(9, •) XC. YC. RC 

I Writ* circulation data for lat*r graphic*. 

DO 20 I “ 0, NPANEL 
WRITEO0.*) I. GS(I). C(I). OG(I) 

20 CONTINUE 

WRITE(10,*) GC, RC 

800 RETURN 

END 


SUBROUTINE RUNCiEI 

Thla program la a mod If lad Runga-Kutto intagratlon Bahama 

COMMON /PANEL/ X(0:505), Y(0:303>, DX(0:S03), DY(0 : S05), 

I W(0;503), DG(0;303) 

COI.WON /TIME/ NPANEU, NSTEP, NTIME, TIME, TSTEP 
COMMON /CONY/ *J(0:505) , V(0:S05) 

COMMON /CORE/ XC. YC, CC, RQ, UC. VC 

/J^/ TFAC, GFAC I TFAC - FACTOR FOR SIZE OF TIME STEP 

WMWN /ZOT/ ZOT 


I Tomporory Baffara 

DIMENSION X0(0:S05), Y0(0:305). DX0(0i305), DY0(0:505), 
# U0(0:505). V0(0:505) 


CALL VELOCT 
TSTEP - , 0025/TFAC 


20 


DO 20 I - 0. NPANEL 

x«(i) - x(i) 

YO(I) - Y(I) 

DX0(I) - OX ( I ) 

DY0(I) - DY(I) 

U0(I) - U(l) 

V0(I) - V(I) 

X(I) - X(I) + U(I) • TSTEP 
Y(I) “ Y(I) + V(I) • TSTEP 
CONTINUE 


1 Vatocltlaa to gat to tha Intarmadlot# point*. 
I Tha Intarmadlot# point* or* ot half-miy. 

I Cat to tha intarmadlot* point*. 


25 


xce - XC 
YC0 - YC 

uce - uc 

VC0 - VC 

XC ■ XC + UC • TSTEP 
YC » YC + VC • TSTEP 

DO 25 I - t , NPANEL-1 

DX(I) - X<I+1) - X( I— 1 ) 

DY(I) - Y(l+t ) - Y(I-1) 

W(I) - SORT ((X( 1+1 ) - X(I-1))**2 + (Y(I+1 ) 
CONTINUE 


i Eatimat* intarmadlot# panala. 


T(I-1))*.2)/4. 


I - NPANEL 


I Ext rape lata ot tha tip. 


(( Y(I-2) - Y( 1-1 ) ) • 
- ( Y(I-1) - Y{ I) ) a 
+ ( Y(I-2) - Y(I~1) ) • 


Y(I) 


8 - 


( X(I- 
( X(I- 
( Y(I- 
( Y( I- 
< Y( I* 
( T(I- 


( Y( 1—1 ) 

/(( X( I-"2) - X(I-1 ) 

- ( x(i~i) - X(I) 

(( X(I-1)..2 - X(I)**2 ) + ( 

- 2. a A * ( X(I-1) - X(I) )) 
DX( I) « - ( Y(I) - B ) 

DY( I) - ( X{ I) - A ) 

W(I) - SORT ( (X(I ) - X(l-1))aa2 + 


- X(I)a.2 ) 

•2)**2 - X(I-1)a*2 ) 
“1)**2 - Y(I)a*2 ) 

■2)*.2 - Y(I”1 )a*2 )) 
-I)-Y(I) ) 

•2) - Y(I~1) )) / 4. 
Y(I=1)aa2 - Y(I)a#2 ) 

/ 2 / ( Y(I-1) - Y(I) ) 


(Y(I) - Y(I~1))a#2)/2, 


CALL VELOCT 

00 35 I • 1, NPANEL 

U(D - ( U0(I) ♦ U(l) ) / 2. 


I Tha mld-cou rat corraction 
I Avarog# In tha mld-coura# corraction 
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V(I) - ( V0<I) + V(I) ) / 2. 

35 continue: 

ue m ( uce + ec ) / 2. 

VC » ( VCO + VC ) / 2. 

TSTEP - .(BK5/TFAC 

DO 50 I * 1, NPANEL I Final translation. 

X(I) - X0;I) •+ U(I) • TSTEP 
V(I) - Yfe(I) + V(I) * TSTEP 
50 CONTINUE 

XC - XC0 + UC * TSTEP 
YC - YC0 + VC • TSTEP 

I Stretching of tip core: conserve volume of the vortex ring. 

RC ■ RC • XC0/XC 

x(0) • 0 

Y(0) - Y{1) 

V(0) - V(1 ) 

TIME - TIME + TSTEP 


RETURN 

END 
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SUBROUTINE EULER 

Thin program (a a Runga-Kutta Intagratlon achama baaed an limiting 
tha largaat convactlon of panala to aoma fraction of tholr wldthe. 

COMMON /PANEL/ X(0: 305), Y(0:505), DX(0;505), DY(0:305), 

§ W(0:305) , DG(0:805) 

COMMON /TIME/ NPANEL, NSTEP, NTIME, TIME, TSTEP 
COMMON /CON V/ U(0:50S), V(0:505) 

COMMON /TFAC/ TFAC , CFAC 

COMMON /CORE/ XC, YC, OC, RC, UC, VC 

CALL VELOCT 

TSTEP - .005 / TFAC 

_ Ponol width baaed time step option .. . 

TSTEP - . 1 

DO 10 1 » 1, NPANEL 

TRAN - W(I) / SQRT( U(I)*«2 + V(I)*«2) / TFAC 
10 IF (TRAN .LT. TSTEP) TSTEP » TPAN 


DO 20 1 - 0, NPANEL 
X(I) - X(I) + U(I) • TSTEP 
Y(I) - Y(I ) + V(I) * TSTEP 
20 CONTINUE 

XC0 - XC 

XC - XC + UC • TSTEP 
YC - YC + VC • TSTEP 

RC * RC • XC0 / XC 

X(0) - 0 
Y(0) - Y(1) 

V(0) - V(1) 

TIME TIME ♦ TSTEP 

RETURN 

END 


SUBROUTINE VELQCT 
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I Thl* program calculate* th* velocity field 
I Induced by th* pans Ik for the program TOLLCALL.FOR 


COMMON /PANEL/ X(0s50t>), Y(0:505). DX(0:305), DY(0:505), 

| W(0:S05), 00(0:503) 

COMMON /CON V/ U(0:S05), V(0:805) 

COMMON /TIME/ NPANEL, NSTEP. NT1ME, TIME, TSTEP 
COMMON /VELO/ R. XX, YY, OXT, DYT, WW, 0 
COMMON /ZQT/ ZOT 

COMMON /CORE/ XC. YC, CC, RC, UC. VC 
PI * ASIN(1.) • 2. 

lllltnminillllUillM Loop over th* troniloUd part*!* 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
DO 10 J " 1, NPANEL 

UJ - 0.0 
VJ - 0.0 

I- Loop over th* inducing panel* - - — 

DO 20 I - t , NPANEL 
R - X(l) 

C - DC(I) 

WW - W(I) 

DXT » DX(I) 

DYT - DY(I) 

XX • X(J) - X(I) 

YY - V(J) - Y(I) 

I ••• Th* Axleymmetrlc caee 
IF (Z0T.CE.1) THEN 

I Interpolate th* e*l (“Induced velocity 

IF ((XX.EO.0).AND.(YY.EQ.0)) THEN 
XX - ~.t • W1 
CALL VELIN ( UINA, VINA) 

CALL VELIO ( UIOA, VIOA) 

CALL VELOUT ( UOUTA. VOUTA) 

XX - .1 • WN 

CALL VELIN ( UINB, VINB) 

CALL VELIO ( UIOB. VIOB) 

CALL VELOUT ( UOUTB. VOUTB) 

UIN - (UINA 4 UINB) / 2. 

VIN ■ (VINA 4 VINB) / 2. 

UIO « (UIOA 4 UIOB) / 2. 

VIO ■ (VIOA 4 VIOB) / 2. 

UOUT * (UOUTA 4 UOUTB) / 2. 

VOUT - (VOUTA 4 VOUTB) / 2. 

ELSE 

CALL VELIN ( UIN, VIN) 

CALL VELIO ( UIO. VIO) 

CALL VELOUT ( UOUT. VOUT) 

END IF 

UJ - UJ 4 UIN 4 UOUT - UIO 
VJ ■ VJ 4 VIN 4 VOUT - VIO 

I ••• 2D Coe*: includ** th* mirror Image component *** 

ELSE 

CALL VEL2D (U1.V1) 

XX « X(J) 4 X(I) 

YY - Y(J) 4 Y(I) 

C • - D0( I ) 

DYT « - DY(I) 

CALL VEL2D (U2.V2) 


UJ ~ UJ + U1 + U2 



VJ « VJ + VI + V2 



END IF 



20 CONTINUE 



|, — — — — — — End of Loop 

U(J) - UJ 

' ■■--■-I 


V(J) « VJ 



10 CONTINUE 



III 1 1 II! 1 1 1 1 1 1 1 1 III 1 1 il 1 III 1 1 1 End of Loop 1 1 1 1 II 1 1 1 1 II 1 1 III II 1 1 HI 1 1 1 1 1 1 1 

min 


1 Colulote the velocity Induced by tho tip-core 



U2 - 0 



V2 0 



R ■> XC 



DO 100 J - 1, NPANEL 



R - XC 



XX - X(J) - XC 
YY » Y(J) - YC 



r. * C «• CC 



IF (ZOT.GE.1) THEN 1 Axleymmetrlc Coee 

S— ■ ww * .5 



i ** CALL VELOUT (U1 , VI) 



: • CO TO 99 



END IF 



CALL VEL2DF (U1 , VI) 1 2D cot* with mirror image 

| XX - X(J) + XC 



c - ~cc 

• 


; CALL VEL2DF (U2. V2) 



99 U(J) - U(J) + U1 + U2 



V(J) - V(J) + VI + V2 

1 


100 CONTINUE 



1 Calculate the velocity Induced on the tip core by the panele 



UC - 0 



VC - 0 



DO 110 I ■ 1, NPANEL 



XX « XC - X(I) 



: V YY - YC - Y( I ) 

DXT - DX(I) 



DYT - DY(I) 
WW»#(I) 

C •» DG(I) 



R » X(I) 



IF (ZOT.GE.1) THEN 1 Axleymmetric Caee 

CALL VELIN (UIN, VIN) 



CALL VELIO (UIO, VIO) 



| CALL VELOUT(UOUT,VOUT) 



U1 - WIN + UOUT - UIO 



VI - VIN + VOUT - VIO 



U2 - 0 



V2 » 0 



GO TO 105 





100 


END IP 

GALL VLL2D (U1.V1) I 2D coaa with mirror Imago 

XX - XC + X(I) 

DYT « DY{I) 

WH - - #(I) 

C « -DG(I) 

CAUL VEL2D (U2,V2) 

IBS UC « UC + 01 + U2 

VC - VC + VI + V2 
110 CONTINUE 

I 2D Com Only: Velocity Induced by tha mirror Imago of tha tip cora 
IF (ZOT.EQ 0) THEN 
XX - 2*XC 
YY - 0.0 
C « ~GC 

CALL VEL2DF (U0.V0) 

UC " UC + U0 
VC * VC + V0 
END IF 


I Axloymmatrlc Coaa Only: Tha Sal f-lnducad Valoclty of tha Cora 
IF (20T.CE.1) VC - VC - CC/4/Pl/XC*( L00(8*XC/RC) - .25) 

U(0) - 0.0 



RETURN 

END 


•a 
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SUBROUT I NE RLO 1 SK ( R , Z . G . N I N , NOIIT , GMAX , SUML , BE LT AL) 

C— THIS DOES CURVE FITTING USING BLENDED PARAUQLA5 ON R, Z, OGAM 
I Madlfled from SPINAL. FOR hy John Kantell* for tie* with panel* 

COMMON /TIME/ NP, NS, NT. T. TS 
COMMON /CO m/ XC, YC, GC 
COMMON /TFAC/ TFAC, Of AC 
COMMON /ZOT/ ZOT 

DIMENSION R(0:305),Z(0:305),G(0:503),S(0i505),QR(0:305), 

$ QZ(0 ; 503) , A0(0: 505) , A1 (0 ; 505) , A2(0 : 505) , A3(r-i 505) , 

$ B0(0: 505) , 01 (0:505) .02(0:505) ,03(0 : 505) ,C0(O : 505) , 

$ Cl (0:505), C2(0:505),C3(0:505),GMID(0:505),5MID(0:505), 

$ 0(0:505) ,RS(0:50D0),ZS(0:5030),AL(0:5050).GS(0:5050) 

DATA NSUB/5/ 

C--GET STRAIGHT LINE DISTANCES BETWEEN POINTS, S(I) IS THE TOTAL 
C DISTANCE ALONG THE STRAIGHT LINES TO NODE I 

S(0>*0.0 
DO 20 I-1.NIN 

20 S(I)-S(I-1)+SQRT((Z(1)»Z(I-1))**2+(R(I)~R(I~1))**2) 


C~ Get value of • for the trailing extended eegrrient, thle value 
C of S l» out at the presumed end of the eLeet. 

S(NIN+1)«S(NIN)+0.3*(S(NIN)~S(NIN-1)) 

C-HJET SLOPES AT EACH INTERIOR NODE, QR(I)-D(R)/DS AT NODE I 

DO 30 I»1,NIN-1 

51- (S(I+1)-S(I))/(S(I»1)”S(1)) 

52- 1 .0/S1 

S3"1 .0/(S(I+1)-S(I-1)) 

QR( I )-( (R( 1-1 )™R( I ) )»S1“(R( 1+1 )-tt( I ) )*S2) «S3 
30 QZ(I)»((Z(I-1)“2(I))*S1-(Z(I+1 )-Z(l))*S2)*S3 

C—COEF'S FOR INTERIOR SEGMENTS 

DO 40 I-1.NIN-2 
S1»1.0/(S(I+1)=S(I)) 

A0( I )-«(!) 

A1(I)-QR(I) 

v A2(I)*S1«("QR(I+1)-2»QR(I)+3»S1*(R(I+1 )«4t(I))) 

A3( I )-S1 *S1 *(QR(I+1 )*0R( I )-2*S1 • (R( 1+1 H»( I ) ) ) 

00(I)»Z(I) 

B1(I)«QZ(I) 

B2( I )*»S1 • ("0Z( 1+1 )-2*QZ( 1 )+3*S1 • (Z( 1+1 )-Z(l ) ) ) 

40 B3(I)"S1*S1*(QZ(I+1)+Q2(I )""2*S1 *(2(1+1 )-Z(I) )) 

C— COEF'S FOR FIRST SEGMENT 

S1»S(1)~S(0) 

S2-S(2)~S(0) 

S3^1 .0/(S(2)-S(l)) 

A0(0)«R(0) 

A1(0)"*((R(1)™R(0))*S2/S1"(R(2)'-R(0))*S1/S2)*S3 

A2(0)-((R(2)"R(0))/S2"(R(1)-R(0))/S1).S3 

A3(0)->0.0 

B0(0)«Z(0) 

B1(0)"((Z(1)-Z(0))*S2/S1~(Z(2)~Z(0))*S1/S2)*S3 


m 


B2(6M(z(2) z(0))/S2-(z(i) Z(e))/si)*r.3 
B3(e)-a.e 

C — COLE'S FOR FAST SEGMENT 


■i i 


N-NIN 

si«s(nHs(n-i) 

S2*S(N-’2) f -S(N-i ) 

S3^1 .0/(S(N-’2)~S(N)) 

A0(fM)Ht(»M) 

A1 (N-1 )«( (R(N)-R(N-I ) ) *52/SMR<N-2M*(N-1 ) ) »S1/S2) *S3 
A2(M-1 )«< <R(N-2)'-RC»+- 1 ))/S2~(R(N)-R(N-’1 ))/S1 )*S3 
A3(N-1)H).0 
B0(N-'1)«Z(fM ) 

B1(N-lH(Z(N)-Z(N-'1))*S2/31-(Z(N-2)-2(rf-l))*S1/82)*S3 
B2(N->1 )-( (Z(N-2)~Z(N^1 ))/S2-(Z(N)~Z(hM ) )/S1 )*S3 
B3(N-1)«0.0 

C — Cst the ongls of th« ring with Indsx N SMALL wrt ths tip cors 
CALL RULANG(NIN.THNIN) 

C™Gst ths "A" cocf f I c I ant for ths Kadln spiral 
RNU-XNU(THNIN) 

IF (A8S(SIN(THNIN) ) .GE. 0.5) THEN 
A»(YC"Z(NIN))/(THNIN**("RNU)*SIN(THNIN)) 

ELSE 

A-(XCHt(NIN))/(THMN*»(^NU)»COS(TMNIN)) 

END IF 

C—GET ARRAY OF R AND Z AT THE SUB;NTERVAL LOCATIONS, AND COMPUTE 
C THE TOTAL LENGTH OF THE CURVE 

C— First do ths blsndsd parabolc ssgmsnts 

RS(0)«A0(0) 

ZS(O)-B0(0) 

AL(0)-0.0 
DO 50 I-0.NIN-1 
DO 50 KSUB-1 .NSUB 
J-I*NSUB+KSUB 

T»(S(I+1)-S(I))*FLOAT(KSUB)/FLOAT(NSUB) 

RS(J)»A0(I)+T*(At(I)+T*(A2(I)+T*A3(I))) 

ZS( J )-B0( I )+T* (B 1 < I )+T* (B2< I )+T«B3{ I ) ) ) 

50 AL(J)-AL(J~1)+SQRT((RS(J)HtS(J-1))**2+(ZS(J)"ZS(J-1))»*2) 

C— Ths nsxt fsw ssctlons find ths volus of thsta at ths snd of ths 
C trolling sxtsndsd ssgmsnt, l.s., gsts ths thsta corresponding 
C to arclsngth 5(NIN+1) 

C — First, bsfors doing ths numsrlcal Intsgrotlon. nssd to gst a 
C volus to uss for dslto thsta 

ITERS-0 

DSGWH(S(NIN+1)-S(NIN))/20.0 

DTG-0.01 

62 TTEMP-THNIN+DTG/2.0 
RNU-XNU(TTEMP) 

ITERS-ITERS+1 

IF (ITERS .CT. 100) STOP ‘ITERS .GT. 160 IN LOOP 61 * 
DSG*A*TTEMF** (-RNtM . 6) *S0RT(RNU**2+TTEMP»*2)«0TG 
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DTGIKPSGW/DGG)*DTU 
If (APS (PI GN- DIG) .LT. 6.061) GO TO B1 
PTG-DTGN 
GO TO 62 
61 PTG-DTGN 

C — Nuw have delta theta (l.e., pit*) which ehould yield gpprnx. 26 
C step* In the numerical Integration from S(N1N) to S(NIN-ll). 

C' — Now do the numerical Integration and ■•urch for theta(NINtl) 

ITERS-0 

SENt>-S(N INtl )-S(NIN) 

STI-H0.0 

TtHTHNIN 

64 RN>XNU(TH) 

ITERS-ITERS+1 

IF (ITEMS .GT. 100) STOP ’ITERS .6T. 100 LOOP 63’ 
ARG1«A»THee(*-RNlM ,0)*SQRT(RNU*«2+TH**2) 

RNU-XNU(TH+DTG) 

ARC2*A*(TH+DTG)»e(>-RNU-1 .0)*SQRT(RNU»«2+(TH»DTG)**2) 

0S»(ARG1+ARG2) *DTG/2 . 0 

IF (STH+DS .GE, SEND) GO TO 63 

STH-STH+DS 

TtNTH+DTC 

GO TO 64 

63 FRAC"(SENB-STH)/DS 
THNIN1-TH+FRAC*DTG 

C'—Now can get the eub-lnlervgl location* over the extended spiral 
C eegmeut. 

DTM THN IN1-THN IN)/NSUB 
00 65 KSUB»1 .NSUB 
THETA*THN I N+KSUB*DTH 
J«NIN*NSU@H<SUB 
RNU-XNU (THETA ) 

RS(J)»XC“A*THETA**(=RNU)*COS(TNETA) 

2S(J)-YC”A«THETA**(-RNU)*SIN(THETA) 

AL(d)-AL(J-1)4SQRT((RS(J)»RS(J-l)).*2+(ZS(J)”ZS(J-1))..2) 

65 CONTINUE 

I Redistribute the Circulation over the Sub-Segments 

DO 751 I • 0, NIN 
DO 751 J • 1, NSUB 
K " I • NSUB + J 
IF (I.EQ.0) THEN 

IF (ZQT.LE.1) THEN 

CS(K) - FLOAT( J )/NSUB • (1.4G(1))/2. - 1. 

GS(0) « -1 . 

ELSE 

GS(K) - FLOAT (J)/NSUB * C(1)/2. 

GS(0) - 0 
END IF 

ELSE IF (I.EQ.NIN) THEN 

DGN - (C(NIN) - G(NIN-1))/2 

GS(K) » FLOAT (J)/NSUB * DGM 4 G(NIN-I) + DGN 

ELSE 

DGI » (0(1+1) - G(I~1))/2 
DGIM1 - (C(I) - 0(1-1 ))/2 

GS(K) ~ (FLOAT(J)/NSUB • DGI) + G(I-1> + DGIM1 


:t? 
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END IF 

continue: 


j. 

V V 
;] 



Q — TOTAE EtNGTH OF CURVE IS SI*;. 

SUME-At«NINF1)*NSUf1) 

1 1 1 M 1 1 1 1 1 1 1 1 1 M M I • 1 1 1 H 1 • M H 11 1 H 1 1 1 1 1 1 1 1 M »• *1 1 1 II 11 1 1 1 1 H H 1 1 1 1 1 1 1 1 II 1 1 • I * 
1 P«n*l Dlncr«tlxatlon Crll»r I hb 

UNO - (NINFl) • NSU0 I Umltn Hnlotlvw OIbIooc* End/Cor# 

XCSEU « SQRT((XC-R5(IEND))**2+(VC-ZS(IEND))**2)/GF , AG 
IF (XESEG.GT.DELTAL) 

| XCSEG « DEITAE I Alwayn Shrinking CrIUrla 

1 1 1 1 1 1 1 1 1 1 II 1 1 II II I U 1 1 H 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 * • •• • * • • 1 * * * 1 * * * * * 1 • * 1 ****** 1 1 * * * 1 * * 
C-- GET NEW R AND Z POINTS 


310 1-1 


00 


DO 60 J«1 ,NSUB*(NINt1) 

Coleulat* thw Pittance to thw po*ltlon« of th« Ruo-punals 

~ 1 ^ _ v a. _ A w > A « \ N / I A \ \ I 


t WM t I W BMW p i - • • • r 

WJM1 - SQRT(( RS(«M)»R(I-l) )*»2 + ( ZS(J-1)~Z(|-1) )**2 ) 
WJM0 * SQRT(( RS(J) - R(*~1) )**2 + ( ZS(J) - Z(l”1) )*»2 > 


) • 


I 


IF (XLSEG .ET. WJM1 .OR. XESEG .GT. WJM0) GO TO 60 
FRAO( XESEG - WJM1 )/( WJM0 - WdMI 
R(I )«FRAO (RS(d)-RS(J-l ))+RS( J-1 ) 

Z( I )-FRAC.*(ZS( J)-ZS(J-1 ) )+ZS{J»1 ) 
G(I)-PRAG*(CS(J)”CS(J~1))#6S(J~1) 

I-I+1 

IF (I.GT.300) GO TO 300 
GO TO 00 
60 CONTINUE 
70 CONTINUE 

R(l) - RS(NSUB*(NIN+1» 

Z{I) - (NSUB»(NSNF1 )) 

G(I) «■ GS(NSUB*(NIM-1)) 

nout - i 

RETURN 


I If fl.>t onoiigF 


Colculota whart 
pon«l odflij will 


ttv 




fiBW 

jp. 
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SUBROUTINE PULANG(NIN.THNIN) 
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I Special Subrout ln« for rout In* REOISK which •xtrapolat'a 
I th* ahaot tdga ualng Kadan'a axponantlal spiral nodal 

COMMON /CORE/ XC, YC, CC 

CQhNQN /PANEL/ R(6:305) .Z<0:505) ,DX(0:505) ,OY(0:309) , 
$ W<0:305). GAMMA(0:505) 

DIMENSION BETA(9:505) 

REAL*4 MAG 
C 

RT-XC 

ZT-YC 

C 

A— 1.0 

B-0.0 

OR(0)~RT 

D*2(0)-ZT 

DQT-A»C+B*D 

CRQSS-A*D~fi»C 

MAG-SORT ( (A»At8*B) * (C»C+0*D) ) 

ARG-OOT/MAG 

IF (ABS(ARG) ,GT. 1.6) ARG»SIGN(1 .O.ARG) 

BLTA(0)-SIGN(1 , 0, CROSS) »ACOS(ARG) 

C 

DO 100 I-1.NIN 
A-R(I-1 )~RT 
8-Z(I=1 )~ZT 
C-R(I)-RT 
D-2( I )"2T 
D0T-A»C+8*D 
CR0SS-A*D~6*C 

MAG-SQRT{ (A*A+B*B)*(C*C+0*D) ) 

ARG-DQT/MAG 

IF (ABS(ARG) .CT. 1.0) ARG-I » IGN( 1 . 0 , ARG ) 

100 BETA(I)-8ETA(I— 1 )+SIGN(1 ,0,CROSS)*ACOS(ARG) 

TMNIN-BETA(NIN) 

RETURN 

END 

FUNCTION XNU(THETA) 

DATA PI/3.141592654/ 

XNU-2. 0/3.0 

IF (THETA .LT. 2.0«PI) THEN 
XNl>-(1.0/(2.e*PI»1 .0))~(1.0/(THETA~1.0))+2. 0/3.0 
END IF 

RETURN 

F.ND 


ktKS 


SUBROUTINE PANELER (X.Y.DX.DY.W.DG.O.NPANEL) 

I Tltl* program takaa th* roult* from routine REDISK and Croat** pan*!* 

DIMENSION X(0:*3S), Y(0:303), DX(0:805), DY(®:303), 

# W(0:503), 00(0:303), 

# 0(0:303), 03(0:303) 

DIMENSION R(0:305), 2(0:303) 

DO 10 I » 1, NPANEL 

R(D - ( X(!) + X(I~1) ) / 2. 

2(1) - ( Y(I) + Y(I~1) ) / 2. 

DX(I) » X(I) - X(I~1) 

DY(I) - Y(I) - Y(I~1) 

W(I) - SQRT< DX(I)*«2 + DY(I)*»2)/2. 

0G(O - ( C(I) - C(M)) / W(I) / 2. 

10 CONTINUE 

DO 20 I - 1 , NPANEL 
X(I) - R(I) 

20 Y(I) » Z(I) 

RETURN 

END 


* 
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Appendix D 

Inputs Used to Generate the Figures in the Text 


Figure 5-7: Elliptical ly Loaded Wing Wake Case 

NPANEL - 80 
TFAC - 1 
GFAC - 1 
NCORE - 4 
20T - 0 

Figure 5-13: Formulation of Vortex Ring by a Disk 

NPANEL - 80 
TFAC - 1 
GFAC - 1 
NCORE - 1 
ZQT - 1 

> 

Figure 5-21: Modelled Helicopter Rotor Wake Case 

NPANEL - 80 
TFAC - 1 
GFAC » 1 
NCORE - 1 
ZOT - 2 

Figure 5-28: Demonstration of Kelvin-Helmlioltz-type Instability 

This experiment -was performed using an early version 
of the simulation program with different set of inputs. 

Figure 5-29: Instability in the Roll-Up of Disk Wake 

NPANEL - 80 
TFAC - 1 
GFAC - 2 
NCORE - 1 
ZOT - 1 
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